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SUMMARY 


A  method  is  presented  for  calibration  of  strain-gauge  balances 
which  does  not  require  that  the  components  can  be  loaded  independently. 
Applicable  to  both  ‘internal’  and  ‘external’  types  of  balance,  the  procedure 
uses  a  single  varying  calibration  load  to  determine  all  linear  and  non-linear 
calibration  coefficients.  Constant  ‘secondary’  loads  on  one  or  more  compo¬ 
nents  are  unnecessary,  although  they  may  be  used  if  desired. 

The  usual  iterative  solution  of  the  second  order  balance  equations 
is  outlined,  and  an  approximate  non-iterative  scheme  is  included  for  com¬ 
pleteness,  though  not  recommended.  Two  methods  of  accounting  for 
dependency  of  the  calibration  coefficients  on  the  signs  of  the  component 
loads  are  presented. 

A  concept  of  ‘buoyancy’  is  introduced  to  simplify  the  application 
of  force  balance  tares,  and  a  procedure  for  determining  the  component  out¬ 
puts  for  absolute  zero  load  (the  ‘buoyant’  offsets)  is  given.  Balance  data  at  a 
series  of  model  attitudes  are  used  to  define  these  offsets,  and  also  the  co¬ 
efficients  in  the  equations  defining  the  component  load  distribution  of  the 
tare  weight  at  any  attitude. 


The  topics  covered  are  ideally  suited  to  formulation  and  solution 
by  matrix  methods,  which  have  been  used  throughout. 


RESUME 


II  s’agit  d’une  methode  d’etalonnage  de  balances  extensometriques 
qui  n’exige  pas  le  chargement  separe  des  composantes.  Applicable  aux 
balances  ‘internes*  et  ‘extemes’,  la  methode  utilise  une  seule  charge  d’etalon¬ 
nage  variable  pour  determiner  tous  les  coefficients  d’etalonnage,  lineaires  et 
quadratiques.  II  n’est  pas  necessaire  d’appliquer  des  charges  ‘secondaries’ 
constantes  sur  une  ou  plusieurs  des  composantes,  mais  on  peut  les  employer 
si  on  veut. 


La  solution  iterative  courante  des  equations  quadratiques  de  la 
balance  est  decrite,  tout  comme  une  autre  methode  approximative,  non 
iterative,  qui  n’est  toutefois  pas  recommandee.  Deux  methodes  visant  a 
compenser  l’effet  du  signe  des  charges  des  composantes  sur  les  coefficients 
d’etalonnage  sont  presentees. 

Le  concept  de  ‘poussee’  permet  de  simplifier  1 ’application  des  tares 
d’equilibrage,  et  une  methode  pour  determiner  la  sortie  des  composantes 
correspondant  a  une  charge  totale  nulle  (corrections  de  ‘poussee’)  est  donnee. 
Des  donnees  experimentales  pour  differentes  orientations  de  la  balance 
servent  a  determiner  ces  corrections  ainsi  que  les  coefficients  des  equations 
de  repartition  de  la  tare  dans  les  composantes,  pour  une  orientation  quel- 
conque  de  la  balance. 

Les  sujets  abordes  sont  particulierement  bien  adaptes,  quant  a 
leur  formulation  et  a  leur  solution,  aux  methodes  matricielles  employees 
tout  au  cours  de  la  presente  etude. 
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A  COMPARISON  OF  METHODS  FOR  CALIBRATION  AND  USE  OF 
MULTI-COMPONENT  STRAIN  GAUGE  WIND  TUNNEL  BALANCES 


1.0  INTRODUCTION 

Extensive  amounts  of  wind  tunnel  testing  time  are  devoted  to  measurement  of  the  static 
aerodynamic  forces  and  moments  acting  on  the  test  model1 .  Careful  calibration  of  the  balances  used  to 
measure  these  loads  is  necessary  to  obtain  the  best  possible  accuracy,  the  objective  being  the  determi¬ 
nation  of  the  constants  in  the  equations  chosen  to  represent  the  actual  behaviour  of  the  balance.  In  a 
test  environment  these  equations  are  solved  to  provide  the  forces  and  moments  corresponding  to  the 
recorded  measurements.  Methods  for  balance  calibration  and  data  reduction  are  discussed  in  this 
report. 


Depending  on  their  location,  wind  tunnel  balances  can  be  classified  into  two  types:  ‘external’ 
balances  which  are  located  outside  the  model  and  test  section,  and  ‘internal’  balances  which  are  located 
inside  the  model  or  its  supports,  or  may  be  integral  with  the  model  (or  a  portion  of  it)  or  the  support. 
The  ‘external’  type  is  generally  regarded  as  being  synonymous  with  the  ‘mechanical’  balances  with 
which  low  speed  wind  tunnels  are  usually  equipped.  In  these,  the  total  aerodynamic  force  and  moment 
are  resolved  into  components  with  the  aid  of  kinematic  mechanisms  consisting  of  links  which  are 
considered  undeformable2,  and  it  is  possible  to  eliminate  almost  completely  interactions  between  the 
components.  However,  as  the  balance  is  outside  the  test  section  the  supports  connecting  the  model  to 
it  are  acted  upon  by  aerodynamic  forces  and  these  effects,  together  with  the  interaction  (interference) 
between  the  model  and  the  supports,  have  to  be  taken  into  account  when  determining  the  true  aero¬ 
dynamic  forces  and  moments  acting  on  the  model.  The  difficulty  of  correcting  for  these  effects  when 
the  magnitudes  increase  at  higher  test  velocities  makes  the  ‘internal’  type  of  balance  more  suitable  for 
high  speed  wind  tunnel  testing. 

‘Internal’  balances  enable  the  influence  of  the  support  to  be  greatly  reduced.  When  the 
balance  is  located  inside  the  model  itself  the  loads  acting  on  the  support  are  not  measured,  although 
the  interference  effect  of  the  support  is  still  reflected  in  the  measured  model  forces.  Such  balances 
find  greatest  application  in  testing  at  high  subsonic,  transonic,  and  supersonic  velocities.  As  wind 
tunnels  having  operating  envelopes  covering  one  or  more  of  these  regimes  typically  have  relatively 
small  test  sections,  the  models  must  be  small  also  and  ‘internal’  balances  are  consequently  synonymous 
with  strain  gauge  measurement  techniques.  In  such  balances  the  kinematic  hinges  of  the  ‘mechanical’ 
type  are  replaced  by  elastic  hinges,  thus  converting  the  ‘metric’  (model)  side  of  the  balance  into  a  kind 
of  floating  frame  connected  to  the  ground  side,  (i.e.  the  model  support),  by  a  statically  determined 
system  of  links.  By  measuring  the  reactions  in  these  links  with  the  aid  of  strain-gauged  elastic  measur¬ 
ing  elements,  the  components  of  the  load  applied  to  the  system  can  be  determined  as  functions  of  the 
strains  in  one  or  more  elastic  elements.  By  suitable  design  of  the  elastic  elements  the  strain  caused  by 
the  component  of  load  or  moment  to  be  measured  can  be  made  to  be  very  much  greater  than  that 
caused  by  any  other  component  of  load.  However,  as  opposed  to  ‘mechanical’  balances,  with  those  of 
the  strain  gauge  type  it  is  usually  not  possible  to  eliminate  completely  the  interactions  between  the 
various  components,  and  these  interactions  must  be  determined  by  the  calibration  procedure. 

In  general  the  process  of  calibrating  a  balance  can  be  defined  as  the  acquisition  of  data 
relating  the  applied  load  to  the  strain  gauge  output  from  each  component,  for  a  sufficient  number  of 
independent  loading  situations  to  allow  determination  of  all  the  constants  in  the  equations  chosen  to 
define  the  behaviour  of  the  balance.  A  common  technique  for  the  calibration  of  ‘internal’  strain  gauge 
balances  involves  the  application  of  loads  which  are  reacted  by  a  single  component  only3,4.  However, 
such  a  procedure  is  not  directly  applicable  to  the  ‘external’  type  where  it  is  either  impossible  to  load 
the  components  independently,  or  to  do  so  would  represent  a  totally  unrealistic  loading  situation  in 
practice.  The  initial  part  of  this  report  describes  a  calibration  technique  which  has  general  applica¬ 
bility  to  both  ‘internal’  and  ‘external’  balances,  and  compares  it  with  the  more  common  technique 
mentioned. 
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The  second  part  of  the  report  is  concerned  with  the  determination  of  the  aerodynamic 
loading  experienced  by  a  model  in  a  wind  tunnel  test.  To  correctly  compute  these  loads  involves, 
firstly,  solution  of  the  balance  equations  to  yield  the  gross  loads  on  the  model,  and  secondly,  calcula¬ 
tion  and  subtraction  of  the  model  weight  tares.  The  common  iterative  approach  to  solution  of  the 
non-linear  balance  equations  is  described,  and  compared  with  an  approximate  method  which  does  not 
require  iteration.  Also,  as  balance  calibrations  often  treat  positive  and  negative  load  behaviour 
separately,  two  methods  of  accounting  for  the  dependency  of  the  calibration  and  interaction  coeffi¬ 
cients  upon  the  sign  of  the  load  are  discussed. 

The  tare  loads,  in  a  balance  axis  system,  are  a  function  of  both  model  weight  and  orientation. 
However,  as  a  result  of  support  deflections,  the  orientation  is  generally  a  function  of  the  gross  load, 
(unless  it  is  sensed  within  the  model  itself),  and  determination  of  the  aerodynamic  loading  normally 
involves  iteration.  A  technique  for  the  application  of  force  balance  tares,  which  avoids  the  iterative 
nature  of  the  problem  by  introducing  a  concept  of  ‘buoyant'  tares,  (tare  readings  for  zero  absolute  load 
on  each  component),  is  described  herein.  The  concept  of  ‘buoyant*  tares  has  the  additional  advantage 
of  avoiding  the  problem  of  ‘initial  load  effects’  for  systems  in  which  the  coefficients  of  second  or 
higher  order  terms  are  non-zero,  as  the  same  origin,  (i.e.  zero  absolute  load),  can  be  used  for  both 
calibration  and  use  of  the  balance. 


2.0  STRAIN  GAUGE  BALANCE  CALIBRATION 

While  a  primary  objective  in  the  design  of  strain  gauge  balances  is  the  minimization  of  inter¬ 
actions  between  components,  i.e.  output  on  a  component  resulting  from  load  applied  to  a  different 
component,  it  is  not  generally  possible  to  eliminate  them  completely.  The  interactions  may  be 
classified  as  either  ‘linear1  or  ‘non-linear’  according  to  whether  they  are  functions  of  a  single  load 
component  or  a  combination  of  load  components.  Typically  the  linear  terms  result  from  such  things 
as  construction  errors  (manufacturing  tolerances  and  assembly  misalignments  for  multi-part  balances), 
improper  positioning  of  the  strain  gauges,  variations  of  gauge  factor,  and  the  electrical  circuits;  the 
non-linear  terms  are  attributable  to  deflections.  In  cases  where  the  design  loads  are  relatively  large  in 
relation  to  the  size  of  the  balance  the  non-linear  effects  will  be  more  significant,  and  it  may  be 
necessary  to  account  for  them  in  order  to  achieve  an  acceptable  measurement  accuracy.  The  decision 
of  whether  to  calibrate  a  balance  as  a  linear  or  non-linear  system  should  be  determined  by  whether  or 
not  the  measurement  accuracy  required  in  the  given  application  can  be  achieved  by  ignoring  non-linear 
effects.  If  a  linear  calibration  suffices  there  is  a  considerable  saving  in  the  effort  required  to  perform 
the  calibration,  and  also  some  computational  simplification,  although  this  latter  point  is  relatively 
insignificant  with  modern  computers. 

For  generality  in  describing  procedures  for  the  calibration  of  multi-component  strain  gauge 
balances,  we  will  consider  one  designed  to  measure  ‘n’  components  of  load  and  moment,  and  will 
assume  that  the  outputs  of  the  strain  gauge  bridge  for  each  component  are,  as  a  result  of  interactions, 
functions  of  all  component  loads. 

From  Reference  3,  the  reading  (strain  gauge  bridge  output)  for  the  ‘i’th.  balance  component 
can  be  expressed  as  a  polynomial  function  of  the  form: 


Ri  *  Ci.l  F1  +  Ci.2  F2  + . +  Ci.n  Fn 

+  Cj  n  Fj2  +  Cj  22  F22  + . +  Ci>nn  Fn2 


+  Ci.l2  F\  F2  +  Ci.I3  F1  F3  +  "  "  *  +  Ci.(n-I)n  F( n - 1  >  Fn 

+  ^  1 1 1  Fj3  +  Cj  222  f23  * - +  Cj  1 23  Ej  F2  F3  +  — ad.  inf.  (1) 

The  output  of  course  depends  upon  the  bridge  excitation  voltage,  but  it  is  usual  to  express  all  quanti¬ 
ties  in  terms  of  unit  voltage. 
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There  will  be  ‘n’  equations  of  this  type  (i  =  1  to  n)  in  which  the  various  constants  are  the 
required  calibration  coefficients.  Practical  experience  with  typical  balances  by  a  wide  range  of  users 
has  indicated  that  terms  of  third  and  higher  degree  can  be  ignored  safely,  and  that  second  degree  terms 
are  nearly  always  small  compared  to  the  direct  linear  term  in  each  equation.  In  what  follows  it  will  be 
assumed  therefore  that  each  component  output  is  represented  to  sufficient  accuracy  by  a  second  order 
polynomial  of  the  form: 


Ri  ”  Ci.l  F1  +  Ci.2  F2  + . +  Ci.n  Fn 

+  Ci.ii  Fl2  +  Ci.22F22  + . +  Cj  nn  Fn2 

+  ^i.12  F1  F2  +  +  ^i.(n-l)n  F(n-1)  Fn  (2) 

For  each  component,  the  coefficients  in  Equation  (2)  can  be  classified  as  follows: 

(a)  ‘linear’,  e.g.  C,  j  for  j  =  1,  n 

(b)  ‘load  squared’,  e.g.  C{ for  j  =  1,  n 

(c)  ‘load  cross  product’,  e.g.  C,-  jk  for  j  =  1,  n  and  k  =  (j+1),  n 

The  number  of  coefficients  represented  by  each  classification  is  seen  to  be  [n],  [n],  and  [n(n  l)/2] 
respectively,  resulting  in  a  total  number  of  coefficients  for  each  component  of  [n(n+3)/2].  Thus,  for  a 
second  order  representation  of  a  balance  designed  to  measure  six  components  of  load  and  moment,  the 
equation  for  each  component  strain  gauge  bridge  reading  contains  27  constant  coefficients  which  must 
be  determined  by  calibration.  Frequently  however,  the  coefficients  in  Equation  (2)  are  assumed  to 
depend  upon  the  signs  of  the  various  component  loads,  to  account  for  any  discontinuity  in  the  varia¬ 
tion  of  the  output  with  load  occurring  at  zero  load.  This  doubles  the  number  of  ‘linear’  and  ‘load 
squared’  coefficients  and  quadruples  the  number  of  second  order  ‘load  cross  product’  coefficients;  the 
total  number  of  possible  coefficients  for  each  component  then  becomes  (2n(n+l)]  for  an  ‘n’  component 
balance  —  84  for  a  6-component  balance,  of  which  27  define  the  component  output  for  a  given  distri¬ 
bution  of  signs  among  the  six  components. 

The  objective  in  calibrating  the  balance  is  the  determination  of  as  many  of  the  sensitivity 
and  interaction  coefficients  as  are  necessary  to  achieve  the  required  measurement  accuracy. 


2.1  Calibration  by  Independent  Loading  of  Components 

The  most  commonly  used  technique  for  calibration3,4  of  internal  strain  gauge  balances  relies 
upon  the  capability  of  applying  load  to  each  component  independently.  Thus,  when  a  single  compo¬ 
nent  is  loaded  the  output  from  that  component  is,  by  definition,  free  of  any  interaction  from  the 
remaining  component  loads  which  are  all  zero,  and  the  outputs  from  the  unloaded  components  are, 
again  by  definition,  the  interactions  caused  by  the  load  applied  to  the  single  loaded  component.  The 
determination  of  the  ‘linear’  and  ‘load  squared’  calibration  coefficients  by  this  method  is  thus  very 
straightforward. 

The  determination  of  the  ‘load  cross  product*  coefficients  in  a  second  order  calibration 
requires  additional  ‘double  component’  loadings  where,  with  the  procedure  commonly  used,  the  load 
on  one  component  is  held  constant  at  some  finite  value  while  that  on  the  other  component  is  varied. 
Data  are  acquired  for  different  values  of  the  constant  load,  from  which  the  ‘cross  product*  coefficient 
can  be  obtained. 

Consider  firstly  the  case  where  only  one  balance  component,  Fj  say,  is  loaded.  The  output 
from  the  strain  gauge  bridge  for  the  Tth.  component  will  then  be  given  by: 


(3) 
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Ri 


Ci,  Fj  +  C, 


j-jj 


Fj2 


Differentiation  w.r.t.  the  variable  load  yields: 


dRi/dFj  =  Cj  j  +  2  CUj  Fj 


(4) 


from  which  it  is  seen  that  the  ‘linear’  coefficient  (C,  ■ )  is  obtained  as: 


Cj.j  =  (3Ri/9Fj)p.  =  0 


(5) 


Further  differentiation  gives: 


a2Ri/9Fj2  =  2Cj.jj  (6) 

from  which  the  ‘load  squared’  coefficient  (Cj  jj)  is: 

Cj.jj  =  i(32Rl/3Fj2)  (7) 


Thus,  for  calibrations  involving  the  -application  of  a  single  load  component,  the  ‘linear’  and 
‘load  squared’  coefficients  are  obtained  directly  by  fitting  a  second  order  polynomial  to  the  variation 
of  each  component  output  with  the  applied  load;  the  required  calibration  constants  are  the  coefficients 
of  the  first  and  second  degree  terms  in  the  applied  load  respectively,  as  indicated  by  Equation  (3). 

To  determine  the  interaction  terms  involving  second  order  cross  products  the  two  components 
involved  must  be  loaded  simultaneously.  If  this  is  accomplished  by  maintaining  a  constant  load  on 
component  ‘k’  (Fk  say)  while  varying  the  load  Fj  on  component  ‘j’,  then  the  equation  defining  the  out¬ 
put  of  component  Y  will  be: 

Ri  =  Cj.j  Fj  +  Cj.k  Fk  +  Cj.jj  Fj2  +  Cj.kk  Fk2  +  Cj.jk  Fj  Fk  (8) 

and  differentiation  w.r.t.  the  variable  load  Fj  now  yields: 

dRj/dFj  =  Cj.j  +  2  Cj.jj  Fj  +  CiJk  Fk  (9) 

or 

(9Rj/3Fj)|;  j=0  =  Cj.j  +  Cj  jk  Fk  (10) 

The  presence  of  the  constant  load  Fk  thus  modifies  the  coefficient  of  the  first  degree  term 
in  the  second  order  polynomial  relating  each  component  output  to  the  variable  load.  The  ‘load  cross 
product’  coefficient  is  isolated  by  further  differentiation  w.r.t.  the  load  Fk,  requiring  the  acquisition  of 
data  at  several  different  values  of  the  constant  load  Fk. 
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Hence: 


Ci.jk  =  a'Ri/aFj  9Fk 


(ii) 


Thus,  provided  the  balance  design  permits  the  application  of  calibration  loads  to  each 
component  independently,  the  procedure  outlined  above  will  allow  determination  of  all  constant 
coefficients  in  the  basic  balance  equation  —  Equation  (2). 

2.1.1  Restrictions  of  the  Independent  Loading  Method 

As  indicated  in  the  Introduction,  independent  loading  of  components  is  not  possible  with 
some  balance  designs  which  should  be  classified  as  ‘external’  rather  than  ‘internal’.  Examples  are 
balances  designed  to  measure  loads  on  ‘half-models’  or  ‘2 -dimensional’  models,  where  the  balance 
components  are  located  outside  the  test  section.  In  such  cases,  even  were  it  possible,  to  apply  load  to 
a  single  component  would  obviously  represent  an  impossible  practical  loading  situation,  contravening 
the  basic  principle  that  the  load  distribution  generated  by  calibration  should  be  representative  of  that 
anticipated  in  application. 

With  the  development  of  the  High  Reynolds  Number  Two-Dimensional  Test  Section  for  the 
NAE  5  ft.  X  5  ft.  Wind  Tunnel,  and  the  associated  sidewall  balance  system,  an  alternative  approach 
was  required  for  calibration  of  this  balance  in  which  no  element  could  be  loaded  independently  of  all 
others.  The  approach  developed  is  completely  general,  the  case  of  independently  loaded  components 
described  above  being  a  special  case. 

2.2  A  Generalized  Calibration  Procedure 

The  basic  premise  in  the  procedure  to  be  described  is  that  the  calibration  consists  of  applying 
a  single  varying  load  to  the  balance,  in  contrast  to  the  method  described  above  in  which  simultaneous 
application  of  multiple  loads  was  necessary  to  obtain  the  non-linear  load  cross  product  terms.  The 
location  and  direction  of  this  single  load  defines  the  distribution  of  the  applied  load  among  the  com¬ 
ponents.  Whether  this  total  (resultant)  varying  load  results  in  load  being  applied  to  one  or  all  of  the 
balance  components  is  immaterial,  although  obviously  more  than  one  component  must  be  loaded  if  the 
non-linear  load  cross  product  terms  are  to  be  evaluated.  Provided  data  are  obtained  for  a  sufficient 
number  of  independent  loading  cases,  substitution  of  these  data  into  the  general  balance  equation  will 
enable  solution  for  the  calibration  coefficients. 

The  strain  gauge  bridge  output  for  the  ‘i’th.  component  is  that  given  by  Equation  (2), 
restated  here  for  convenience: 


Ri  =  C, 

.1  F, 

+  Cj.2 

f2  + . 

+  C,„Fn 

+ 

C-i.i  i 

F,2  + 

Ci.22  F22  +  -  ' 

. +  Ci.nn  Fn2 

+ 

Ci.12 

F,  F2 

(12) 

The  output  (Rj)  can  also  be  expressed  as  a  second  order  polynomial  in  the  resultant  varying  applied 
load  (P)  thus: 

Ri  =  roi  +  a(  P  +  bj  P2  (13) 


Differentation  of  (12)  w.r.t.  the  load  P  yields: 
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Thus, 


0F 


l  3F2 


0Rj/9P  =  Cj  .  -  +  C:  2  -  + 

'  >1  gp  >-2  gp 


0Fj 

^  2  0:  II  Fl 

111  1  gp 


9F2 

^  2  Ci  22  F2  ~  ^ 

..22  2  gp 


+  c 


i.(n-l)n 


9F; 


+  Ci.l2 

9F/ 


9F2  0F, 

-  +  F, - 

9P  2  0P 


(n-l) 


n  (n-l) 

-  +  F  - 

9P  "  9P 


9Fj  9F2 

(9Rj/9P)P=0  =  ai  =  Cj.,—  +  Ci.2—  + 


+  C: 


9Fn 

9P 


(14) 


(15) 


and  for  a  fixed  location  of  the  applied  load  we  may  write: 

j  dp  p 


representing  the  fraction  of  the  applied  load  P  which  is  reacted  by  the  ‘j’th.  balance  component.  For  a 
fixed  position  of  the  applied  load  these  ‘component  load  distributions’  or  ‘load  ratios’  are  constant  and 
defined  by  the  static  equilibrium  of  the  system. 

Equation  (15)  therefore  becomes: 


fl  +  Ci  2f2  ♦ . +  ci>nfn 


(16) 


and,  considering  all  balance  component  outputs,  there  are  ‘n’  equations  of  this  type  each  containing 
‘n’  unknown  coefficients.  By  acquiring  data  for  a  minimum  of  ‘n’  independent  loading  cases, 
(locations  of  the  applied  load),  the  resulting  linear  simultaneous  equations  can  be  solved  to  provide 
the  linear  calibration  coefficients. 

To  define  the  non-linear  coefficients,  further  differentiation  of  (14)  w.r.t.  P  gives: 


/9F,\2 

92Ri/9P2  -  2CiI1\ap‘/ 
/9F| 

'U2\9P 


+  2C: 


or 


92Rj/9P2  =  bj  «  Cj. 


11  M 


f.2  +  C. 


i.22  ‘2 


+  Ci.l2  f2  + 


+  Ci.(n-l)n  f(n-l)  fn 


(17) 
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Again,  the  definition  of  a  sufficient  number  of  independent  equations,  [n(n+l)/2]  in  this  case,  by  the 
application  of  the  load  P  at  various  locations  and  orientations,  will  yield  sets  of  linear  simultaneous 
equations  which  can  be  solved  to  provide  the  ‘load  squared’  and  ‘load  cross  product’  calibration 
constants. 

2.2.1  Matrix  Formulation  of  the  Equations 

The  sets  of  linear  simultaneous  equations  defined  by  Equations  (16)  and  (17)  are  most 
conveniently  represented  in  matrix  notation, 

[S]  =  [C]  [L]  (18) 

(n,m)  (n,m)  (m,m) 


where 

n  is  the  number  of  balance  components, 

m  is  the  number  of  calibration  coefficients  per  component,  [m  =  n(n+3)/2] 

[S]  contains  a  combination  of  the  curve  fit  coefficients  for  each  loading  case,  ‘n’  of  the  first 
degree  and  [n(n+l)/2]  of  the  second  degree  in  the  applied  load 

[C]  contains  the  calibration  coefficients,  and 

[L]  defines  the  distribution  of  the  applied  load  P  between  the  balance  components  for  each 
loading  case. 

For  simplicity  in  what  follows,  the  number  of  balance  components  can  be  taken  as  three 
without  loss  of  generality.  In  referring  to  the  curve  fit  coefficients  and  the  component  load  ratios,  the 
first  subscript  will  serve  to  define  the  component  and  the  second  the  data  for  a  particular  loading  case. 
For  the  balance  calibration  coefficients  the  subscripts  refer  to  the  ‘interacted -upon’  and  ‘interacting’ 
components  respectively  in  the  normal  way. 

For  the  case  of  the  six  independent  loading  cases  necessary  to  determine  all  calibration 
coefficients  for  a  three-component  balance,  the  elements  of  the  three  matrices  [S] ,  [C] ,  and  [L]  are: 


[S]  = 

all 

a2 1 

a12 

a22 

al  3 

a23 

bil 

b21 

b12 

b22 

b13 
b2  3 

b14 

b24 

b15 

b25 

b16 

b26 

(3,9) 

_a31 

a32 

a33 

b3i 

b32 

b33 

b34 

b35 

b36_ 

C,.i 

Cl. 2 

C1.3 

Cl. n 

Cl. 22 

C1.33 

Cl. 12 

Cl. 13 

C|  .23 

[C]  = 

c2.i 

C2.2 

C2.3 

C2.1 1 

C2.22 

C2.33 

C2.12 

C2.1  3 

C2.23 

(3,9) 

_p3.1 

C3.2 

C3.3 

C3.11 

C3.22 

C3.33 

C3.12 

C3.13 

C3.23 

_ 

"*11 

*12 

*13 

0 

0 

0 

0 

0 

0 

“ 

[LI  = 

*21 

*3. 

0 

0 

*22 

*32 

0 

0 

*23 

*33 

0 

0 

0 

0 

*.l2 
*21 2 

0 

0 

*122 

*222 

0 

0 

*132 

*232 

0 

0 

*142 

*242 

0 

0 

*152 

*252 

0 

0 

f.62 

*262 

(9,9) 

0 

0 

0 

*3. 2 

*322 

*332 

f  2 
*34 

*35  2 

*36  2 

0 

0 

0 

*11  *21 

*12  *22 

*13  *23 

*14  *24 

*15  *25 

*16  *26 

0 

0 

0 

*1!  *31 

*12  *32 

*13  *33 

*14  *34 

*15  *35 

*16  *36 

0 

0 

0 

*21  *31 

*22  *32 

*23  *33 

*24  *34 

*25  *35 

*26  *36 
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Equation  (18)  can  be  re-written,  partitioning  the  matrices  thus: 


[  a  j  b  ] 
(3,3)  (3,6) 


[Cl  !  C2  ] 
(3,3)  (3,6) 


LI  I  0 
(3,3)  ' 

“0  ,  L2 

I  (6,6) 


where  the  various  sub-matrices  contain: 

a  —  polynomial  linear  coefficients  (a iz) 
b  —  polynomial  second  order  coefficients  (biz) 

[Cl]  —  linear  calibration  coefficients  (Cjj ) 

[C2]  —  non-linear  calibration  coefficients  (Cj  jj  & 

[LI  ]  —  ‘load  ratios’,  i.e.  the  component  loads  expressed  as  ratios  of  the  applied  load 

[L2]  —  ‘squares’  and  ‘cross  products’  of  the  load  ratios 

and  subscript  V  designates  the  particular  loading  case. 

To  solve  for  the  balance  calibration  coefficients,  Equation  (19)  is  re-written  as: 


IC1IC2J  -  [a[bj  (—!•—]  ' 

fir1  !  o  1 

=  [a  j  b] 

L  o  ;l2-1J 


from  which 


[Cl]  =  [a]  [LI]-1 


[C2]  =  [b]  [L2] 


and  clearly  the  solutions  for  the  linear  and  second  order  coefficients  are  independent.  Provided  that 
the  load  ratio  matrices  [LI]  and  [L2]  are  non-singular,  the  calibration  matrices  are  easily  determined. 

2.2.2  Components  Loaded  Independently 

It  is  readily  shown  that  a  calibration  performed  by  loading  the  balance  components  inde¬ 
pendently  and  singly  is  a  particular  case  of  this  general  solution.  For  this  case  we  have: 


Fii/P  -  f„  -  1 


Fij/P  "  rij  *  0 


Again  using  the  three-component  system  for  ease  of  illustration,  we  have: 


[Cl]  =  [a] 


1  0 
0  1 
0  0 


0 

0 

1 


fa] 


(24) 


indicating  that  the  linear  calibration  coefficients  are  simply  the  coefficients  of  the  first  degree  terms  of 
the  polynomial  relating  component  output  and  applied  load. 

For  the  non-linear  coefficients  only  the  ‘load  squared’  terms  can  be  derived  from  single 
component  loadings,  and  the  matrix  [L2]  is  consequently  reduced  in  size  from  (6,6)  to  (3,3).  Thus: 


[C2] 


[b] 


1  0 
0  1 
0  0 


lb] 


(25) 


indicating  that  the  calibration  coefficients  for  the  ‘load  squared’  terms  are  simply  the  coefficients  of 
the  second  degree  terms  of  the  load/output  polynomial  equation. 

Equations  (24)  and  (25)  simply  express  the  relationships: 


C|.j  =  (dRi/3Fj)Fj  =  Q 

C,jj  =  ~  d2R,/d Fj2 


for  ‘i*  and  *j’  varying  from  1  to  the  number  of  components* 


2.2.3  Effect  of  Constant  Load  on  One  or  More  Components 

The  procedure  described  in  Section  2.1  utilized  the  technique  of  applying  a  constant  load  to 
one  component  in  combination  with  a  variable  load  on  another  in  order  to  determine  the  coefficient 
of  the  second  order  load  cross  product  in  question.  However,  the  formulation  just  described  for 
calibration  using  multi-component  loading,  is  based  on  the  premise  that  the  only  load  applied  to  the 
balance  is  a  single  varying  one  and  that  consequently  each  of  the  individual  component  loads  is 
proportional  to  this  load.  This  will  no  longer  be  true  if  a  constant  load  is  applied  to  one  or  more 
components,  in  addition  to  the  variable  load.  The  formulation  is  extended  below  to  allow  the  use  of 
this  calibration  procedure. 

Assume  again  an  V  component  balance  having  a  second  order  load/output  relationship,  and 
let  the  calibration  procedure  consist  of  applying  a  single  resultant  varying  load  (P)  in  combination  with 
one  or  more  constant  component  loads  (Kif  i  <  *  n). 

In  terms  of  the  balance  calibration  coefficients,  the  output  of  the  ‘i’th.  component  is  now: 
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Ri  =  Cj  i  (F,  +  K,)  +  Cj  j  (F2  +  K2)  + . +  Ci  n  (Fn  +  Kn) 

+  Cj.n(F, +K,)2  +  Ci  22(F2+K2)2  + . +  Ci  nn  (Fn  +  Kn)2 

+  Cj  i2  (F,  +  Kj)  (F2  +  K2)  +  CM3  (F,  +K,)(F3  +Kj)  +  -  --- 

+  ^i.(n-l)n  (^(n-l)  +  K(„-i))  (Fn  +  Kn)  (26) 

and  in  terms  of  the  variable  applied  load  (P)  it  is,  as  before: 

Rj  =  r0j  +  a;  P  +  bj  P2 

Differentiating  Equation  (26)  w.r.t.  the  variable  load  (P)  yields: 

dRj/dP  =  Cj.,  f,  +  Cj  2  f2  + . +  Cj  n  fn 

+  2  Oj  |  j  (F,  +  K, )  f,  + . +  2  Cj  nn  (Fn  +  Kn)  fn 

+  Cj  I2  [(F,  +  Kj)  f2  +  (F2  +  K2 )  f i ]  + . 

+  Cj.(n.1)n  [(F(n_i)  +  !<(„_„)  fn  +  (Fn  +  Kn)  t(n.u]  (27) 

and  hence 

3Rj 

ai  =  Tp"  =  Ci.i  fj  +  +  Cj  n  fn  +  2Cji,fIK]  + - 

or  p=o 

+  2  Cj  nn  fn  Kn  +  Cj  j  2  (*i  K2  +f2  K,)  + - 

+  ^i.(n-l)n  (f(n-l)  +  K(n-i))  (28) 

Further  differentiation  w.r.t.  P  gives: 

1  d2Ri 

bi  =  T  —  =  Cj.nf,2  + . +  Cj.nnfn2 

2  0P2 

+  Cj.,2f,f2  +----  +  Cj.(n_1)nf(n.nfn  (29) 


It  can  be  seen  that  the  presence  of  the  constant  loads  affects  only  the  coefficient  of  the  first 
degree  term  in  the  load /output  equation,  and  now  it  is  possible,  in  principle,  to  solve  for  all  calibration 
constants,  linear  and  non-linear,  using  only  the  ‘a’  coefficients  from  the  curve  fitting  of  the  calibration 
data. 
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In  matrix  notation: 


[a]  =  [Cl  !  C2]  [L*]  (30) 

(n,m)  (n^n)  (m,m) 

where  ‘m*  is  the  total  number  of  calibration  coefficients  per  component,  and  the  elements  in  a  column 
of  the  matrix  [L*]  are  functions  of  the  ‘load  ratios’  and  the  constant  loads  as  indicated; 

1^1 »  ^2  *  ^3  *  ‘  *  "  "  ^n» 

2  f  i  K| ,  2  f2  K2,2f3K3,----2fn  Kn, 

[f,  K2  +  f2K,),- ---IV,) Kn  +  fnK<n.1)|} 

Using  the  ‘b’  polynomial  coefficients  we  have,  as  previously, 

[b]  =  [C2]  [L2]  (31) 

(n,p)  (n,p)  (p,p) 


where  ‘p’  is  the  number  of  non-linear  calibration  coefficients  per  component. 

Whilst  in  principle  Equation  (30)  can  be  used  to  determine  all  coefficients,  to  do  so  requires 
the  acquisition  of  data  for  [n(n+3)/2]  independent  loading  situations.  This  requirement  can  be 
reduced  by  [n]  cases  if  Equation  (31)  is  used  to  solve  for  the  second  order  coefficients,  (matrix  [C2] ), 
as  in  this  case  only  [n(n+l)/2]  load  applications  are  required.  By  writing  Equation  (30)  in  the  form: 


*  [Cl)  [LI*]  +  [C2]  [L2*]  (32) 

(n,n)  (n,n)  (n,p)  (p,n) 


where  p  =  (m  -  n), 

we  see  that  by  substituting  the  known  matrix  [C2] ,  the  ‘a’  values  can  be  corrected  for  the  effect  of 
the  constant  loads.  The  linear  calibration  coefficients  are  then  obtainable  from  the  corrected  ‘a’  values 
thus: 

[Cl]  -  |[a]  -  [C2]  [L2*]}  [LI*]'1  (33) 

Solution  of  this  equation  requires  only  V  loading  cases,  which  may  be  chosen  from  the  [n(n+l)/2] 
already  used  to  solve  for  the  second  order  constants. 

If,  as  a  result  of  the  loading  procedure  used,  the  matrix  [L2]  in  Equation  (31)  should  be 
singular,  it  becomes  necessary  to  solve  for  a  matrix  [C2]  containing  only  the  ‘load  squared’  calibration 
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coefficients.  This  situation  will  result  if  the  procedure  described  in  Reference  3  and  Section  2.1  is 
followed,  whereby  the  balance  components  are  loaded  independently.  In  this  case  the  product  of  the 
‘load  ratios'  of  any  two  components  would  always  be  zero,  the  ‘load  ratios’  being  the  ratios  of  the 
single  variable  applied  load  reacted  by  each  component,  and  therefore  unaffected  by  the  presence  of 
constant  loads.  Removal  of  the  load  ratio  product  terms  from  the  matrix  [L2]  would  allow  determi¬ 
nation  of  the  ioad  squared’  calibration  coefficients  by  substitution  of  ‘n’  independent  loading  cases 
into  Equation  (31 ).  The  procedure  for  correcting  the  ‘a’  values  is  still  applicable,  but  the  matrix  [Cl  ] 
in  Equation  (33)  would  now  be  composed  of  both  the  ‘linear’  and  the  ‘load  cross  product’  calibration 
coefficients,  and  [n(n+l)/2]  load  cases  would  be  required  to  achieve  the  solution.  The  matrices  [LI*] 
and  [L2*]  would,  of  course,  be  modified  to  reflect  the  changed  composition  of  [Cl]  and  [C2] . 

2.2.4  Advantages  of  the  Generalized  Procedure 

The  most  obvious  advantage  is  that  it  is  equally  applicable  to  either  ‘internal’  or  ‘external’ 
types  of  balances.  Also  the  formulation  is  such  that  the  calibration  constants  determined  by  the 
procedure  relate  to  the  ‘fundamental’  components  and  not  to  ‘derived’  components.  Here  ‘fundamen¬ 
tal’  components  are  defined  as  those  into  which  the  system  of  elastic  hinges  and  links  resolves  the 
applied  load  and  moment;  ‘derived’  components  are  those  obtainable  from  combinations  of  these 
fundamental  reactions.  For  example,  a  simple  two-component  balance  might  have,  as  fundamental 
components,  the  forces  reacted  at  two  locations  in  the  same  plane.  These  two  forces  also  define  a 
resultant  force  and  a  moment  but,  according  to  the  definition  used  here,  these  are  derived  rather  than 
fundamental  components.  The  importance  of  maintaining  this  distinction  becomes  apparent  when 
considering  the  use  of  separate  calibrations  for  positive  and  negative  loads  in  each  component.  As 
mentioned  before,  this  is  frequently  done  to  compensate  for  discontinuities  in  the  load/output  rela¬ 
tionship  at  zero  load,  however  such  a  distinction  is  meaningful  only  when  considered  in  terms  of  the 
fundamental  components.  For  example,  in  the  case  of  the  two-component  balance  just  cited,  a 
positive  total  resultant  load  (a  derived  component)  indicates  only  that  both  fundamental  component 
loads  cannot  be  negative,  and  likewise  a  particular  sign  of  the  total  moment  may  result  from  three  of 
the  four  possible  sign  combinations  of  the  fundamental  components.  Clearly  the  behaviour  of  the 
balance  depends  on  the  signs  of  the  reactions  in  the  basic  flexure  elements  and  not  on  the  signs  of  the 
total  resultant  force  and  moment. 

The  concept  of  having  multiple  components  loaded  simultaneously  allows  the  calibration 
procedure  to  be  much  more  representative  of  the  manner  in  which  the  balance  will  be  used  in  practice. 
From  a  purely  theoretical  point  of  view,  provided  that  the  actual  behaviour  is  exactly  as  defined  by 
the  chosen  equation  for  the  load/output  relationship,  the  precise  combination  of  loading  arrangements 
used  in  calibration  is  immaterial;  all  would  yield  exactly  the  same  result.  However,  in  practice  the 
chosen  polynomial  relationship,  whether  it  be  limited  to  second  order  or  includes  higher  order  terms, 
will  not  describe  exactly  the  behaviour  of  each  component,  and  this  results  in  the  solution  of  the 
system  of  equations  being  dependent  upon  the  conditions  represented  by  the  equations.  In  other 
words,  the  calibration  coefficients  will  depend  upon  the  manner  in  which  the  balance  was  loaded 
during  calibration,  as  different  loading  distributions  will  result  in  different  deflection  characteristics. 
This  being  so,  there  would  appear  to  be  a  strong  argument  in  favour  of  applying,  during  calibration, 
loads  which  generate  distributions  representative  of  those  likely  to  be  experienced  in  use  of  the  balance 
to  measure  model  forces  in  a  wind  tunnel.  The  occurrence  of  single  component  loadings  in  a  practical 
test  situation  would  be  very  infrequent,  and  accordingly  the  preferred  method  of  calibration  would 
seem  to  be  the  proposed  multi-component  loading  procedure. 

This  procedure  may  also  enable  a  satisfactory  accuracy  to  be  obtained  without  the  need  to 
make  a  complete  non-linear  calibration,  something  which  can  realize  a  considerable  saving  in  calibra¬ 
tion  effort.  This  comes  about  by  virtue  of  the  fact  that,  during  a  calibration  involving  the  loading  of 
several  components  simultaneously,  the  strain  gauge  bridge  outputs  automatically  include  the  effects, 
linear  and  non-linear,  of  all  components  having  non-zero  loads.  If  these  outputs  are  now  ‘linearized’ 
by  straight  line  curve  fitting  of  the  output  data  with  respect  to  the  applied  load,  the  ‘slopes’  will  take 
into  account,  in  a  linear  fashion,  all  of  the  ‘load  squared’  and  'load  cross  product’  non-linear  effects  of 
those  components  with  non-zero  load.  In  the  case  of  a  single  component  loading  only  the  ‘load 
squared’  effect  of  the  loaded  element  would  be  present.  Provided  that  the  non-linear  effects  are 
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reasonably  small,  it  is  possible  that  a  linear  calibration  utilizing  compound  loading  of  the  components 
may  provide  adequate  accuracy,  better  than  that  obtainable  from  a  linear  calibration  determined  from 
single  component  loadings. 

2.2.5  An  Optimized  Solution  of  the  Calibration  Equations 

As  the  calibration  coefficients  will  depend  upon  the  particular  set  of  equations  solved  (for  a 
practical  situation),  and  as  the  formulation  permits  the  use  of  calibration  loadings  chosen  to  be  repre¬ 
sentative  of  the  expected  test  measurement  environment,  this  procedure  lends  itself  to  the  idea  of 
obtaining  an  optimum  solution  from  a  large  number  of  calibration  loadings.  These  could  be  made  to 
cover,  as  completely  as  possible,  the  expected  magnitude  and  location  of  the  test  load,  and  the  result¬ 
ing  overdetermined  system  of  equations  could  then  be  solved,  in  a  least  squares  fashion5 ,  to  provide 
the  optimum  solution. 

Assume  that  we  wish  the  solution  to  an  overdetermined  system  of  equations,  of  the  form  of 
Equation  (30),  for  the  (n,m)  matrix  of  calibration  coefficients  for  an  ‘n’-component  balance,  utilizing 
‘z*  loading  cases  where  V  is  greater  than  ‘m\  and  ‘m’=[n(n+3)/2] . 

We  therefore  require  to  solve  for  the  linear  [Cl]  and  non-linear  [C2]  calibration  matrices 
in: 


[a]  -  [Cl  !  C2]  [L*]  (34) 

(n,z)  (n,m)  (m,z) 


This  overdetermined  system  of  equations  can  be  solved,  in  a  least  squares  sense,  by  post-multiplication 
of  each  side  of  Equation  (34)  by  the  transpose  of  matrix  [L*] .  Thus: 


[a]  [L*]t  =  [Cl  !  C2]  [L*]  [L*]T  (35) 

(n,z)  (z,m)  (n,m)  (m,z)  (z,m) 

and  hence 

[Cl  !  C2]  *  [a]  [L*]T[[L*J  [L*]T1  '  (36) 

(n,m)  (n,z)  (z,m)  [(m,z)  (z,m)  J 

(m^n) 


In  thia  manner  many  different  loading  distributions  may  be  permitted  to  influence  the  resulting  cali¬ 
bration  coefficients,  which  represent  a  least  squares  solution  of  all  the  equations  (load  cases)  involved; 
weighting  factors  may  be  applied  to  particular  loading  configurations  if  desired. 

2.3  Some  Notes  on  ‘Zeros’  and  Intercepts  in  Curve  Fitting 

Fundamental  to  all  methods  of  reducing  calibration  data  to  the  form  of  sensitivity  and 
interaction  constants,  is  the  process  of  obtaining  the  coefficients  in  the  polynomial  load /output 
equation  which  provide  a  best  fit  to  the  experimental  data.  In  the  foregoing,  the  equation  chosen  has 
been  a  second  order  polynomial  for  the  general  case,  although  obviously  a  straight  line  fit  would  be 
used  if  a  linear  calibration  was  the  objective. 

The  presence  of  other  than  first  degree  terms  presents  the  problem  of  defining  what  the  true 
'zero'  on  the  load  scale  should  be  as,  for  anything  other  than  a  linear  relationship,  the  use  of  one  ‘zero’ 
for  calibration  and  another  for  application  will  give  rise  to  errors.  As  indicated  in  the  Introduction, 


. _ tet. 
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this  report  will  describe  a  procedure  for  the  determination  of  ‘buoyant’  (i.e.  zero  absolute  load) 
component  ‘zeros’  for  reduction  of  data  from  wind  tunnel  balances.  Accordingly,  if  the  same  ‘buoyant’ 
condition  is  used  as  the  ‘zero’  during  calibration  of  the  balance  the  potential  problem  of  initial  load 
effects6,7  is  avoided.  With  the  extended  formulation  given  in  Section  2.2.3  this  objective  is  easily 
achieved  by  specifying  the  calibration  tare  loads,  (i.e.  the  weight  of  the  calibration  body  and  the 
‘metric  mass’  of  the  balance  itself),  as  constant  loads  existing  during  the  calibration.  The  precise  dis¬ 
tribution  of  this  tare  weight  between  the  balance  components  will  depend  upon  the  balance  orientation 
and  the  centre  of  gravity  of  the  tare  mass,  and  must  be  determined  for  each  loading  case.  Of  course, 
the  significance  of  the  tare  load  lies  in  its  magnitude  relative  to  the  load  range  over  which  the  balance 
is  to  be  used,  and  also  in  the  relative  importance  of  non-linear  effects  for  the  balance;  a  consideration 
of  these  facts  will  indicate  whether  or  not  the  tare  load  need  be  considered  for  a  particular  circum¬ 
stance. 


The  preceeding  discussion  of  ‘zeros’  was  concerned  with  zeros  on  the  load  scale.  The  zeros 
for  the  strain  gauge  bridge  outputs  are  the  intercepts  (rG)  in  Equation  (13),  and  are  of  no  real  signifi¬ 
cance  as  balance  data  reduction  is  concerned  with  incremental  changes  and  the  intercept  values  cancel 
in  the  process  of  tare  subtraction.  The  main  significance  of  the  intercept  is  that,  in  general,  the  curve 
fitting  should  be  performed  with  a  ‘free*  intercept  in  order  to  obtain  the  best  fit  to  the  experimental 
data.  Having  found  the  coefficients  in  Equation  (13)  —  *r0\  ‘a’,  and  ‘b*  —  which  define  this  best  fit, 
the  intercept  values  may  be  disregarded. 

There  is  one  exception  to  this  general  rule,  namely  the  fitting  of  a  straight  line  to  data  which 
contains  some  non-linear  effects.  An  example  would  be  in  making  a  linear  analysis  of  calibration  data 
from  multi-component  loadings,  as  suggested  in  Section  2.2.4. 


Suppose  that  the  experimental  output  is  well  fitted  by  the  second  order  polynomial: 

R  »  r0,  +  a,  P  +  b,  P2 

as  indicated  in  the  diagrams  above.  If  this  same  data  is  fitted,  using  a  least  squares  technique,  by  the 
straight  line: 


R  ■  ro2  ♦  a2  P 

the  result  will  be  that  the  sum  of  the  deviations  of  the  data  on  either  side  of  the  line  will  cancel.  The 
greatest  deviations  will  therefore  occur  at  the  end  points  and  in  the  middle  of  the  range  as  shown  in 
the  left  hand  diagram. 
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Now  consider  both  of  these  curves  translated  along  the  ‘output'  axis  so  as  to  pass  through 
the  origin,  as  indicated  by  the  dashed  lines  in  the  diagram  —  this  is  the  equivalent  of  subtracting  a  tare 
value.  We  see  that  the  straight  line  now  gives  a  maximum  deviation  from  the  data  at  the  middle  of  the 
range,  and  that  all  points  on  it  are  biased  in  one  direction  relative  to  the  experimental  data.  This  is 
clearly  an  undesirable  situation,  but  it  can  be  corrected  by  forcing  the  linear  fit  to  pass  through  the 
same  origin  as  the  second  order  fit;  i.e.  by  fitting  the  data  to  the  straight  line: 


R  *  rQi  +  a3  P 


In  practice  this  is  accomplished  by  subtracting  the  value  rOJ  from  all  experimental  outputs  and  then 
fitting  with  the  straight  line: 


R  -  &3  P 


.  For  a  least  squares  fit  of  this  equation  we  must  minimize  the  sum  of  the  squares  of  the 
residuals  w.r.t.  the  coefficient  ‘aj’;  i.e. 


2  [Rj-  a3  Pj]2 


Differentiating  w.r.t.  ^3’  and  equating  to  zero: 


2  2  Pj  [Rj  -  a3  Pi]  =  0 

indicating  that  the  deviations  from  the  straight  line,  when  weighted  by  Pj,  will  sum  to  zero;  this  will 
bias  the  line  toward  data  at  large  values  of  P4  which  is  undesirable  —  line  [1]  in  the  right  hand  diagram. 

To  obtain  a  better  fit,  in  which  the  deviations  will  simply  sum  to  zero,  it  is  necessary  to 
weight  each  point  by  a  factor  inversely  proportional  to  the  square  root  of  the  load  ?r  Thus  we  require 
to  minimize: 


r  1  v 

—  <R.  -  *3  Pi) 

p.m 


Differentiating  w.r.t.  ‘a3’  and  equating  to  zero: 


2  2  [Rj  -  a3  Pj]  -  0 


and  the  required  slope  of  the  straight  line  is  seen  to  be  simply: 


a3  *  2  Rj/2  Pj 


This  straight  line  is  such  that  the  deviations  of  the  data  on  either  side  of  the  line  cancel  —  a  common 
interpretation  of  ‘best  fit’.  A  comparison  of  this,  (line  [2] ),  and  the  original  data  is  given  in  the  right 
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hand  diagram;  it  will  be  noted  that  the  uni-directional  bias  no  longer  exists,  and  that  the  deviation 
generally  increases  as  the  load  increases.  This  trend  is  quite  acceptable  as  greater  absolute  errors  can  be 
tolerated  when  the  load  is  large. 

Following  a  recent  discussion*  on  this  topic  it  was  suggested  to  the  author  that  a  logical 
extension  to  the  approach  of  balancing  the  deviations  on  either  side  of  the  fitted  straight  line  would  be 
to  balance  the  deviations  when  expressed  as  fractions  of  the  local  data  value.  This  would  come  closest 
to  the  ‘ideal’  of  error  being  a  constant  percentage  of  measurement. 

To  achieve  such  a  ‘fit’  each  point  must  be  weighted  by  the  factor: 

l/(Rj  Pi)I/2 


and  we  require  to  minimize 


(Rs  P>),/2 


(Ri-  «3  Pi) 


Differentiating  w.r.t.  ‘a3’  and  equating  to  zero: 


22 


-  0 


from  which  the  coefficient  a3  is  seen  to  be: 


a3  -  N/2 


for  N  data  values.  The  straight  line  for  which  the  slope  (a3)  is  defined  in  this  way,  (line  [3] ),  would 
seem  to  be  preferable  to  the  one  for  which  the  unweighted  deviations  sum  to  zero. 

Thus,  contrary  to  what  might  be  anticipated,  in  making  a  linear  curve  fit  of  experimental 
balance  calibration  data,  the  best  technique  is  to  force  the  straight  line  to  pass  through  the  experimen¬ 
tal  output  value  at  zero  load,  or  alternatively,  through  the  intercept  value  obtained  from  a  second 
order  curve  fit  of  the  data.  This  will  eliminate  the  uni-directional  bias  resulting  from  the  use  of  a  ‘free’ 
intercept  in  the  curve  fit,  combined  with  the  subtraction  of  ‘tares’  in  reducing  the  data  which  results 
in  cancellation  of  the  common  intercept  value.  In  addition,  suitable  weighting  of  the  data  can  be  used 
to  achieve  a  balance  of  the  deviations,  expressed  as  percentage  of  measurement,  on  either  side  of  the 
line. 


*  with  S.A.  Sjolander  of  DSMA  International  Inc.,  Toronto. 
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2.4  Normalization  of  the  Calibration  Matrix 

As  the  constants  in  the  linear  [Cl]  and  non-linear  [C2]  calibration  matrices  are  closely 
related  to  the  polynomial  coefficients  of  the  first  and  second  degree  in  the  load  respectively,  it  is 
obvious  that  they  are  not  non-dimensional.  The  linear  constants  have  units  of  ‘output  unit  per  unit  of 
load’,  while  those  for  the  non-linear  constants  are  ‘output  unit  per  unit  of  load  squared’.  Consequently 
the  complete  matrix  of  calibration  constants  is  applicable  only  when  the  balance  is  used  with  a  system 
having  the  same  output  units  as  the  calibration  system. 

Rather  than  having  to  re-specify  all  of  the  constants  to  suit  a  different  output  system,  a 
‘normalized’  form  of  the  matrix  can  be  used,  in  which  the  constants  in  the  equation  for  each  compo¬ 
nent  are  expressed  as  ratios  of  the  primary  sensitivity  of  the  respective  component.  These  primary 
sensitivities  are  the  diagonal  elements  of  the  linear  calibration  matrix  [Cl] ,  and,  following  the 
normalization  procedure,  are  the  only  elements  of  the  calibration  matrix  which  retain  a  dependence 
on  the  system  used  to  measure  the  strain  gauge  bridge  outputs.  A  factor  for  converting  from  one 
output  system  to  another  can  be  obtained  readily  by  either  check  loading  individual  components,  (so 
as  to  exclude  interactions),  or  by  unbalancing  the  individual  component  strain  gauge  bridges  using 
shunt  resistors.  The  primary  sensitivities  are  multiplied  (or  divided),  by  this  factor  to  convert  them 
from  one  system  to  the  other. 

Normalization  is  achieved  by  division  of  the  constants  in  the  equation  for  each  component 
by  the  primary  sensitivity  of  that  component.  This  corresponds  to  pre-multiplication  of  the  original 
calibration  matrix  by  a  diagonal  matrix  [D]  composed  of  the  reciprocals  of  the  diagonal  elements  of 
the  linear  matrix  [Cl] .  The  resulting  ‘interaction’  matrix  [X]  is  then: 


[X]  =  [XI  !  X2]  =  [D]  [Cl  !  C2]  (37) 

(n,m)  (n,n)  (n,p)  (n,n)  (n,n)  (n,p) 


where  the  elements  of  [X]  have  universal  application,  but  those  of  [D]  are  dependent  on  the  particular 
readout  system. 

The  ultimate  check  of  the  complete  calibration  procedure  can  now  be  made  by  applying 
check  loads  to  the  balance  and  using  the  derived  calibration  coefficients  to  re-calculate  the  loads  from 
the  measured  component  outputs.  Several  procedures  for  solution  of  the  balance  equations  will  now 
be  compared. 


3.0  REDUCTION  OF  STRAIN  GAUGE  BALANCE  DATA 

When  the  balance  components  are  loaded,  each  component  strain  gauge  bridge  exhibits  a 
reading,  related  to  all  of  the  components  acting,  which  is  defined  by  the  equation: 


R  =  [Cl]  F  +  [C2]  F* 


(38) 


where 

R  is  the  vector  of  outputs  at  unit  excitation  voltage, 

[Cl]  is  the  matrix  of  linear  calibration  coefficients, 

[C2]  is  the  matrix  of  non-linear  calibration  coefficients, 

F  is  the  vector  of  component  loads, 

and 

F*  is  a  vector  composed  of  the  squares  and  cross-products  of  the  component  loads. 
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A  common  procedure  for  determination  of  the  force  vector  F  in  this  equation  is  given  in 
Reference  3.  This  form  of  solution,  which  requires  iteration,  is  restated  here  as  it  is  the  correct 
procedure  fundamentally.  It  is  compared  with  an  approximate  solution  procedure  known  to  have 
been  used,  and  perhaps  still  in  use,  at  some  wind  tunnel  test  facilities  in  Europe.  This  latter  approach 
has  in  its  favour  the  fact  that  the  approximation  assumption  avoids  the  need  for  an  iterative  type  of 
solution;  whether  this  simplification  justifies  the  approximation  involved  is  another  question. 

Pre-multiplication  of  both  sides  of  Equation  (38)  by  the  diagonal  matrix  [D]  in  Equation  (37) 

gives: 


[D]  R  =  [D]  [Cl]  F  +  [D]  [C2]  F* 


or 


F~  =  [XI]  F  +  [X2]  F* 


(39) 


where  the  vector  F~  consists  of  the  approximate  component  loads  (ignoring  interactions)  resulting 
from  division  of  each  component  output  by  the  appropriate  primary  sensitivity  constant. 

Solving  for  the  force  vector  F  gives: 


F  =  [XI]-'  l~  -  [XI]-1  [X2]  F* 


(40) 


For  balances  having  an  assumed  linear  behaviour  the  solution  is  given  explicitly  by: 

F  =  [XI]-1  F~  (41) 


For  non-linear  behaviour  the  iterative  solution  is  commenced  using  the  solution  of  Equation  (41),  Fj 
say,  to  calculate  a  first  approximation  of  the  non-linear  term  (Aj )  thus: 

A,  =  [XI]-1  [X2]  F,*  (42) 

where  the  elements  of  the  vector  Fj  *  are  composed  of  the  squares  and  cross-products  of  the  elements 
of  the  approximate  vector  Fj . 

A  second  approximation  (F2)  is  then  obtained  as: 


F2  =  F,  -  A,  (43) 

from  which  F2*  may  be  determined  and  then  substituted  into  Equation  (42)  to  provide  a  further 
estimate  of  the  non-linear  term. 

In  general  the  final  solution  is  given  by: 


En  =  El  —  <n- 1 ) 


(44) 


where  ‘n’  is  the  number  of  iterations. 
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Any  iterative  solution  raises  questions  as  to  whether  or  not  successive  approximations  will 
always  converge,  and,  if  they  do,  whether  the  limit  is  a  unique  solution  of  the  given  equations.  These 
questions  are  considered  in  Reference  6  where  it  is  shown  that  a  converged  solution  represents  a 
unique  solution,  and  that  the  question  of  whether  or  not  convergence  takes  place  is  dependent  upon 
the  balance  interaction  coefficients  and  the  design  loads  only.  A  criterion  is  given  which  can  be  evalu¬ 
ated  for  a  given  balance  prior  to  its  use;  for  typical  balances  the  non-linear  effects  are  generally  quite 
small  and  convergence  is  usually  rapid. 

3.1  An  ‘Explicit’  Approximate  Solution 

In  Equation  (40)  the  elements  of  the  non-linear  force  vector  F*  are  functions  of  the  squares 
and  cross-products  of  the  true  load  vector  F,  i.e. 


F*  =  f(F) 


In  the  method  to  be  outlined  it  is  assumed  that  the  non-linear  terms  are  sufficiently  small,  in  compari¬ 
son  with  the  linear  terms,  that  the  non-linear  vector  F*  may  be  approximated  using  the  elements  of 
the  vector  Fj  obtained  by  ignoring  non-linear  effects.  Thus 


F*  ~  F~*  =  f(Fj) 

where 

Fj  =  [XI]-1  F~  (45) 


Conceding  for  the  moment  the  validity  of  this  approximation,  it  is  now  possible  to  substitute 
for  F*,  (in  terms  of  F~ ),  in  Equation  (40),  thus  making  the  solution  for  the  true  load  vector  (F)  explicit 
in  terms  of  F^  and  the  non-linear  vector  F^  *  where: 


F~ *  *  f(F~)  and  F~  =  [D]  R 


as  seen  from  Equation  (39). 

References  8  and  9  describe  analytically  ways  in  which  this  substitution  for  the  non-linear 
vector  can  be  made.  However,  the  following  computational  method  is  offered  here  as  it  provides  a 
means  of  ‘tailoring’  the  assumption  to  particular  load  ranges  and  distributions. 

It  is  required  to  determine  the  coefficient  matrices  in  the  expression: 

F  =  [XI]'1  F~  -  [B]  F~*  (46) 

Restating  Equation  (40)  we  have: 

F  =  [Xl]”1  F~  -  [XI]'1  [X2]  F* 
and  from  a  comparison  of  this  with  Equation  (46)  it  is  clear  that: 


[B]  F~*  =  [XI]- 1  [X2]  F* 


(47) 
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Now  the  interaction  matrices  [XI]  and  [X2]  are  known,  and  for  any  specified  force  vector 
F  we  may  compute  the  output  vector  R  from  Equation  (38)  and  hence  F~ .  Then: 


F*  =  f(F) 

and 

F~*  =  f(F~ ) 

Thus,  by  computing  the  vectors  F*  and  F~  *  corresponding  to  a  sufficient  number  of  input 
load  vectors  F,  the  required  matrix  [B]  is  obtained  from  the  equation: 

[B]  =  [XI]'1  [X2]  [F*]  (48) 

(n,p)  (n,n)  (n,p)  (p,p)  (p,p) 


where  p  =  [n(n+l)/2]  and  ‘n’  is  the  number  of  components. 

The  ‘p’  columns  of  [F*]  and  [F~  *]  are  composed  of  the  corresponding  vectors,  (F*  and  F~  *),  for 
each  input  load  vector  F.  Thus,  by  appropriate  choice  of  these  load  vectors,  the  solution  for  [B]  can 
be  ‘tailored’  to  specific  combinations  and  ranges  of  the  individual  component  loads.  The  value  of  [B] 
determined  from  the  analytic  treatments  of  References  8  and  9  is  obtained  when  the  magnitudes  of 
the  component  loads  in  the  chosen  vectors  are  made  to  approach  zero.  As  ‘p’  input  load  vectors  are 
used  in  the  solution  for  [B] ,  it  follows  that  these  particular  vectors,  i.e.  component  load  distribution 
and  magnitudes,  must  be  exact  solutions  of  Equation  (47),  and  thus  of  both  the  approximate  and  the 
iterative  solution  methods.  Use  of  the  approximate  method  to  compute  load  conditions  close  to  any 
of  these  particular  load  vectors  will  thus  involve  little  error. 

Use  of  ‘q’  input  load  vectors,  (q  >  p),  in  solving  for  [B]  results  in  an  overdetermined  set  of 
equations  which  may  be  solved,  in  a  least  squares  sense5,  by: 


thus  allowing  a  greater  number  of  load  vectors  to  influence  the  result,  although  none  will  be  an  exact 
solution  of  Equation  (47). 

This  approach  to  balance  data  reduction,  in  addition  to  involving  an  approximation,  suffers 
from  the  disadvantage  that  it  is  not  amenable  to  the  use  of  different  positive  load  and  negative  load 
sensitivity  and  interaction  coefficients;  what  might  be  termed  ‘plus-minus’  coefficients.  This  is  easily 
seen  from  the  fact  that  the  product  of  [X2]  and  the  inverse  of  [XI]  is  involved  in  the  determination 
of  [B] ,  which  therefore  will  be  valid  only  for  that  particular  ‘sign  of  load’  distribution.  The  only 
advantage  of  the  approach  seems  to  be  that  no  iteration  is  involved  in  the  solution;  given  the  rapid 
convergence  of  the  iterative  solution  for  typical  balances,  and  the  speed  of  modern  computers,  this 
does  not  seem  to  the  author  to  be  sufficient  justification  for  an  ‘unnecessary’  approximation.  The 
method  is  not  employed  at  NAE  for  reduction  of  balance  data  from  the  5  ft.  X  5  ft.  Wind  Tunnel. 

3.2  Accounting  for  ‘Plus-Minus’  Calibration  Coefficients 

Calibrations  of  strain  gauge  balances  which  separate  the  behaviour  of  the  fundamental 
components,  (i.e.  the  flexure  elements),  to  positive  and  negative  loading  are  quite  common.  The 
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inherent  ‘cubic’  type  of  load/output  relationship  for  typical  bending  beam  flexure  elements  is  fitted 
better  by  a  second  order  polynomial  if  the  positive  and  negative  load  data  are  assessed  separately.  Two 
methods  are  presented  to  account  for  this  dependence  of  the  calibration  coefficients  on  the  sign  of  the 
loads. 

3.2.1  Synthesis  of  a  Specific  Matrix  for  Each  Load  Vector 

The  calibration  coefficients  defining  the  balance  behaviour,  in  number  and  type,  are: 

Tinear’  [2n],  +  and  -  loads 

ioad  squared’  [2n],  +  and  -  loads 

Toad  cross-product’  [2n(n-l)],  ++,  +  -,  -+,  and  --  combinations  of  loads. 

Thus,  in  total,  [2n(n+l)]  coefficients  must  be  stored  for  use  by  the  data  reduction  programme. 

In  solving  for  the  force  vector  corresponding  to  a  given  output  vector  only  [n]  ‘linear’,  [n] 
Toad  squared’,  and  [n(n  l)/2]  Toad  cross-product’  coefficients  are  involved  for  any  given  sign  combi¬ 
nation.  Thus,  if  the  combination  of  signs  in  the  required  load  vector  can  be  established,  then  the 
appropriate  linear  [XI]  and  non-linear  [X2]  matrices  can  be  synthesized  from  the  stored  coefficients, 
and  Equation  (40)  solved  to  yield  the  load  vector.  Note  that  the  solution  requires  the  inversion  of 
matrix  [XI]  and  determination  of  the  product  of  this  inverse  and  the  non-linear  matrix  [X2] . 

The  difficulty  lies  in  establishing  the  correct  sign  distribution  upon  which  to  base  the  choice 
of  the  appropriate  coefficients.  Strictly  these  must  be  chosen  according  to  the  signs  of  the  loads  but 
these  are,  as  yet,  unknown.  The  only  option  is  to  select  the  matrix  elements  on  the  basis  of  the  signs 
of  the  component  outputs,  following  subtraction  of  tares  to  remove  any  ‘offsets’.  However,  depending 
on  the  exact  definition  of  the  tares  used,  this  distribution  of  the  signs  of  the  outputs  may  not  corre¬ 
spond  to  the  distribution  of  the  signs  of  the  loads.  For  example,  subtraction  of  a  tare  corresponding 
to  the  model  weight  will  give  the  indication  of  a  zero  output  when  actually  the  components  which 
react  the  weight  are  negatively  loaded.  It  is  necessary  therefore,  to  compare  the  distribution  of  signs 
in  the  calculated  load  vector  with  that  used  in  selecting  the  matrix  elements;  if  differences  exist  a  new 
selection  of  the  matrix  elements  must  be  made,  using  the  updated  sign  distribution,  and  the  load  vector 
re-computed.  This  procedure  should  be  repeated  until  the  correct  load  distribution  is  found. 

The  use  of  ‘buoyant’  tares,  to  be  described  later,  facilitates  the  determination  of  the  correct 
sign  distribution  by  defining  the  tares  as  the  component  outputs  at  zero  absolute  ioad.  Consequently, 
following  subtraction  of  these  ‘buoyant’  tares,  the  sign  distribution  of  the  output  vector  and  that  of 
the  load  vector  are  generally  identical  and  no  iteration  is  required.  It  is  possible  that  interactions  could 
cause  a  difference  between  the  two  distributions  for  components  which  are  very  lightly  loaded;  how¬ 
ever,  this  situation  would  likely  result  in  failure  of  the  iterative  scheme  to  converge  to  a  unique 
solution,  providing  instead  a  ‘flip-flop’  between  two  load  distributions.  Under  these  circumstances  the 
choice  would  be  arbitrary  in  any  event. 

3.2.2  Concept  of  Load  and  Absolute  Value  of  Load 

In  this  approach  Equation  (39)  is  rewritten  to  express  the  elements  of  the  vectors  F  and  F* 
in  terms  of  the  ‘signed*  and  ‘absolute’  values  of  the  individual  component  loads.  The  linear  term  in 
Equation  (39)  thus  becomes: 


[XI]  F  =  [XL1J  F  +  [XL2]  FL2 


where 


[XL1]  and  [XL2]  are  (n,n)  matrices  appropriate  to  the  ‘signed’  and  ‘absolute  value’  load 
vectors  F  and  FL2  respectively. 
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The  second  order  term  can  be  divided  into  two  ‘load  squared*  and  four  ‘load  cross- 
product*  contributions  thus: 

[X2]  F*  =  [XSQ1]  FSQ1  +  [XSQ2]  FSQ2 
+  [XCP1]  FCP1  +  [XCP2]  FCP2 
+  [XCP3]  FCP3  +  [XCP4]  FCP4 

where 

[XSQ1]  and  [XSQ2]  are  (n,n)  matrices  appropriate  to  the  ‘load  squared*  vectors  FSQ1 
and  FSQ2  respectively,  “ 

[XCP1],  [XCP2],  [XCP3],  and  [XCP4]  are  (n,p)  matrices  appropriate  to  the  ‘load  cross- 
product’  vectors  FCP1,  FCP2,  FCP3,  and  FCP4  respectively, 

and 

p  -  [n(n-l)/2]  as  usual. 

The  general  elements  for  all  of  the  load  vectors  noted  are  as  follows: 

F  :  {(Fj),  i  =  1,  n} 

FL2  :  { (IFj|),  i  =  1,  n } 

FSQ1  :  {(Fj  •  Fj),  i  =  1,  n} 

FSQ2  :  {(Fj  •  |F,|),  i  =  1,  n} 

FCP1  :  |  {(Fj  •  Fj),  i  =  1,  n],  j  =  (i+1),  n  } 

FCP2  :  { [(Fj  •  |Fj  |),  i  =  1,  n] ,  j  =  (i+1),  n } 

FCP3  :  { [(|Fj|  •  Fj),  i  =  1,  n],  j  =  (i+1),  n} 

FCP4  :  { [(|Fj  •  Fj|),  i  =  1,  n],  j  =  (i+1),  n} 

where 

Fj  and  Fj  are  the  component  loads. 

Using  these  vectors  and  calibration  matrices,  Equation  (39)  becomes: 

F~  =  [XL1  ]  F  +  [XL2]  FL2  +  [XSQ1]  FSQ1  +  [XSQ2]  FSQ2 

+  [XCP1]  FCP1  +  [XCP2]  FCP2  +  [XCP3]  FCP3  +  [XCP4]  FCP4  (50) 

and  the  solution  for  the  load  vector  F  is  now  given  by: 

F  =  [XL1]  '  F~  -  [XU]'1  |  [XL2]  FL2  +  [XSQ1]  FSQ1  +  [XSQ2]  FSQ2 

+  [XCP1]  FCP1  +  [XCP2]  FCP2  +  [XCP3]  FCP3  +  [XCP4]  FCP4 )  (51) 
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This  appears  more  complex  than  Equation  (40)  but  it  should  be  noted  that  the  inverse  of 
the  matrix  [XL1] ,  and  the  product  of  this  and  the  other  matrices,  can  be  pre-stored  as  they  are 
independent  of  the  sign  distribution  of  the  load  vector. 

Defining  the  general  elements  of  the  ‘linear',  ‘load  squared',  and  ‘load  cross-product’  coeffi¬ 
cient  matrices  as: 


‘linear'  —  j+,  %l  }~ 

‘load  squared’  —  qi  j+,  qt j" 

‘load  cross-product’  —  q j++,  q  j+",  q  j~+,  qj“~ 

the  expressions  for  the  general  elements  of  each  of  the  derived  matrices  are: 


[XLl] 

:  (V  +  Ci.f )/2 

[XL2] 

••  <V  - 

[XSQ1] 

:  (qM+  +  qi.j")/2 

[XSQ2] 

:  (q,/  -  qj.p/2 

[XCP1] 

:  (Ci.j++  +  c^-  + 

Ci/'  +  ci.j  +)/4 

[XCP2] 

:  <ci.j++  -  c,f  " 

Ci.j+-  ♦  ci.j  +)/4 

[XCP3] 

:  (Ci/+  "  Si"  + 

Cj/-  -  Cj.j -+)/4 

[XCP4] 

••  (Ci.j++  ♦  c,"  - 

Ci.i+“  -  ci.j”+)/4 

3.2.3  Comparison  of  the  Two  Approaches 

Firstly,  for  both  approaches  the  storage  requirement  for  the  coefficient  matrices  is  identical  — 
[2n(n+l)]  coefficient:  per  component  for  a  complete  second  order  system. 

The  first  method  requires,  for  each  load  vector  to  be  determined,  synthesis  of  the  calibration 
matrix  for  the  appropriate  distribution  of  load  signs,  and  a  matrix  inversion  and  multiplication  to 
determine: 

[XI]-1  and  [XI  r1  [X2] 


in  Equation  (40).  This  procedure  may  have  to  be  repeated  in  order  to  arrive  at  the  correct  sign 
distribution  for  synthesis  of  the  matrix,  unless  the  ‘buoyant’  tares  are  utilized.  Solution  for  the  load 
vector  is  explicit  for  a  linear  calibration,  and  iterative  for  a  non-linear  calibration. 

The  second  method  utilizes  eight  pre-stored  matrices,  (two  ‘linear’,  two  ‘load  squared',  and 
four  'load  cross-product’),  to  eliminate  the  requirement  to  select  and  invert  a  matrix  for  each  load 
vector.  However,  a  wider  range  of  load-related  vectors  must  be  calculated,  and  a  larger  number  of 
matrix  multiplications  are  required.  The  solution  is  iterative  for  both  linear  and  non-linear  calibra¬ 
tions. 


T 
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Some  tests  were  made  to  determine  the  relative  computation  time  requirements  of  the  two 
methods,  by  solving  for  a  series  of  load  vectors  (for  a  six -component  system)  using: 

(a)  a  linear  calibration, 

(b)  a  non-linear  calibration,  ignoring  all  second  order  ‘load  cross-product’  terms,  and 

(c)  a  complete  second  order  non-linear  calibration. 

In  all  cases  the  method  involving  synthesis  of  a  new  matrix  for  each  load  vector  proved  to  be  less 
efficient;  for  cases  (a)  and  (b)  the  ratio  of  the  computation  times  was  approximately  3.5,  which  was 
reduced  to  approximately  2  by  inclusion  of  the  ‘load  cross-product’  terms  in  case  (c).  Thus,  if  compu¬ 
tation  time  is  considered  to  be  a  significant  factor,  the  method  involving  the  use  of  separate  pre-stored 
coefficient  matrices  would  seem  to  be  preferable. 


4.0  THE  APPLICATION  OF  FORCE  BALANCE  TARES 

In  a  wind  tunnel  test  measuring  model  forces,  meaningful  aerodynamic  data  is  obtained  only 
after  determination  and  removal  of  the  model  weight  effects  from  the  measured  data.  The  distribution 
of  the  model  weight  tares  among  the  balance  components  is  a  function  of  the  magnitude  and  centre  of 
gravity  position  of  the  ‘metric  mass’  and  the  balance  orientation  relative  to  a  gravitational  axis  system. 
However,  unless  the  orientation  is  sensed  within  the  model/balance  assembly,  as  a  result  of  support 
deflections  it  will  depend  upon  the  total  load  acting.  Thus  the  tare  values  cannot  be  determined  with¬ 
out  knowing  the  orientation,  which  cannot  be  determined  without  knowing  the  total  load,  which,  in 
turn,  cannot  be  determined  without  knowing  the  tare  values,  .  .  . .  ;  the  solution  thus  requires  iteration. 

This  requirement  for  iteration  arises  from  the  fact  that  the  test  (wind-on)  measurements  are 
regarded  as  being  incremental  from  some  finite-load  tare  condition,  for  which  the  component  loads,  in 
a  balance-axes  system,  depend  upon  the  orientation.  If  instead  one  supposes  an  ‘absolute  zero  load’  or 
‘buoyant’  condition  for  the  tare,  (which,  by  definition,  will  be  independent  of  the  orientation,  ‘metric 
mass’  and  centre  of  gravity  position),  then  all  incremental  measurements  relative  to  this  tare  are  pro¬ 
portional  to  load  on  an  absolute  scale.  The  total  load,  (and  hence  the  support  deflection  and  balance 
orientation),  are  thus  available  without  iteration,  and  the  component  load  distribution  can  be  deter¬ 
mined  from  a  knowledge  of  the  magnitude  and  centre  of  gravity  location  of  the  ‘metric  mass’.  The 
concept  of ‘buoyant’  tares  thus  considerably  simplifies  the  determination  of  the  aerodynamic  loading  — 
the  objective  of  the  test. 

The  component  strain  gauge  bridge  ‘buoyant’  offsets  are  hypothetical  of  course  as  the 
balance  is  never  ‘weightless’.  However,  it  will  be  shown  that  they  can  be  derived  from  tare  data  taken 
at  different  balance  orientations;  these  data  also  serve  to  define  the  magnitude  and  centre  of  gravity 
location  of  the  ‘metric  mass’.  The  formulation  given  is  general  and  self-contained;  it  calculates  all 
necessary  information  from  the  tare  (wind-off)  data  obtained  at  a  series  of  balance  orientations,  and 
hence  it  is  unnecessary  to  determine  the  model  weight  and  centre  of  gravity  by  separate  means. 

Another  virtue  of  the  ‘buoyant’  concept  has  been  mentioned  previously  in  connection  with 
calibration  of  the  balance;  namely  the  avoidance  of  ‘initial  load  effects’  for  non-linear  balance  systems 
by  the  use  of  ‘buoyant’  tares  for  both  calibration  and  application  of  the  balance. 

4.1  Balance- Axes  Load  Components  of  the  Metric  Mass 

In  the  concept  of  ‘buoyant’  tares  it  is  necessary  to  consider  the  ‘metric  mass’  as  opposed  to 
the  weight  of  the  model.  This  ‘metric  mass’  is  defined  to  include  all  parts  of  the  model/balance  system 
whose  weight  results  in  loads  on  any  of  the  balance  components,  the  magnitude  and  distribution  of 
which  depend  upon  the  orientation  of  the  balance  relative  to  a  gravitational  system  of  axes.  Due  to 
the  method  of  construction  of  the  balance,  it  is  possible  that  different  components  may  sense  differ¬ 
ent  ‘metric  mass’  contributions  from  parts  of  the  balance,  and  therefore  different  total  tare  weights 
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and  centre  of  gravity  positions.  This  will  be  illustrated  later  by  the  fact  that,  in  a  six -component 
balance  for  example,  expressions  for  the  weight  and  centre  of  gravity  position  are  defined  separately  by 
the  balance  components  measuring  loads  in  the  normal,  side,  and  axial  force  planes.  In  what  follows, 
‘tare  weight’  is  used  to  indicate  the  product  of  the  ‘metric  mass’  and  the  gravitational  constant. 

A  very  common  type  of  model  support  mechanism  utilizes  freedom  in  pitch  (a)  and  roll 
(0)  to  position  the  model  at  any  desired  orientation  in  the  wind  tunnel.  The  components  of  the  tare 
weight,  in  a  system  of  balance  axes,  are  given  in  terms  of  these  angles  by: 


Wx  =  W  sin  a 
Wy  -  W  cos  a  sin  0 
Wz  =  W  cos  cc  cos  0 


where  Wx,  Wy,  and  W2  are  the  tare  weight  components  along  the  ‘x\  *y\  and  V  axes  respectively  —  see 
Figure  1. 


Note  that  the  choice  of  pitch  and  roll  angles  for  definition  of  model  orientation  is  not 
restrictive,  as  the  matrix  of  direction  cosines  defining  the  position  of  the  axes  system  can  be  generated 
by  a  sequence  of  Euler  rotations.  If  the  first  rotation  applied  is  one  about  the  vertical  axis  it  will  have 
no  effect  on  the  component  load  distribution  of  the  tare  weight.  Further  rotations  of  pitch  and  roll 
respectively  will  thus  provide  the  equivalent  pitch  and  roll  angles  for  the  given  balance  orientation, 
regardless  of  the  sequence  of  rotations  used  to  define  the  orientation  originally.  The  determination  of 
different  Euler  rotation  sequences,  corresponding  to  a  given  matrix  of  direction  cosines,  is  described  in 
Reference  10. 

The  net  tare  loads  sensed  by  the  balance  components  can  be  represented  by  the  vector  FN 
which  is  defined  as: 


FN  =  W  [G1  sin  a  +  G2  cos  a  sin  0  +  G3  cos  a  cos  0] 


(52) 


where  the  elements  of  the  vectors  Gl,  G2,  and  G3  are  functions  of  the  geometric  position  of  the 
balance  components  and  the  position  of  the  centre  of  gravity  of  the  tare  load  W.  By  combining  these 
three  vectors  and  the  weight  we  define  a  matrix  [G]  thus: 


|G]  =  IWG1  !  WG2  1  WG3] 


and  can  now  re-write  Equation  (52)  as: 


sin  a 

FN  =  [G]  cos  a  sin  0 
cos  a  cos  0 


(53) 


Extending  this  same  notation  to  express  the  component  loads  for  several  different  tare  orienta¬ 
tions  we  have: 


(54) 
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sin  aj  i 

sin  Oi2  j  . 

. .  .  j  sin  am 

[FN] 

=  [G] 

cosa1sin01  ] 

cos  a2  sin  02  !  • 

•••  [  cos  am  sin  0m 

cos  cos  (j>i  \ 

cos  a2  cos  (f>2  !  . 

.  .  .  |  cos  am  cos  <t>w 

(n,m) 

(n,3) 

(3,m) 

where  ‘n’  is  the  number  of  balance  components,  and 
‘m’  is  the  number  of  tare  attitudes. 

For  convenience  we  can  denote  the  matrix  comprising  the  trigonometric  functions  of  the 
tare  angles  by  [T]  so  that  Equation  (54)  becomes: 


[Fn]  =  [G]  [T]  (55) 

(n,m)  (n,3)  (3,m) 


This  equation  defines  the  balance  component  loads  in  terms  of  the  unknown  matrix  [G] 
and  the  tare  attitude  expressed  by  the  pitch  and  roll  angles.  Obviously  if  the  elements  of  [G]  can  be 
determined  from  the  tare  data,  this  relation  provides  the  means  to  determine  tare  load  effects  at  any 
attitude. 

4.2  Expressions  for  the  Elements  of  [G] 

For  a  particular  design  of  balance  the  elements  of  the  matrix  [G]  can  be  defined  analytically 
from  a  consideration  of  the  static  equilibrium  of  the  system.  For  the  case  of  the  generally  familiar  type 
of  six -component  balance  manufactured  by  the  Able  (formerly  Task)  Corporation,  the  component  load 
vector  is: 


Fn  =  {N1,N2,  Yl,  Y2,  X,  r} 


representing  the  forward  and  aft  normal  and  side  force  components,  and  the  axial  force  and  rolling 
moment  respectively.  The  matrix  [G]  is  given  by: 


[G]  =  W 
(6,3) 


-z0/2a 

z0/2a 

y0/2b 

-y0/2b 

1 

0 


0 

0 

0.5  (1  +  x0/b) 
0.5(1  -  x0/b) 
0 


-0.5(1  +  x0/a) 
-0.5(1  -  xQ/a) 
0 
0 
0 

y0 


where  (x„,  yQ,  z0)  define  the  centre  of  gravity  of  the  tare  weight,  relative  to  the  balance  reference  point 
midway  between  the  forward  and  aft  components, 

2a  is  the  normal  force  component  spacing, 

2b  is  the  side  force  component  spacing, 

and  the  sign  convention  for  x0,  y0,  and  z0  is  positive  forward,  to  starboard,  and  downward 
respectively  —  see  Figure  1. 
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Inspection  of  these  expressions  for  the  elements  of  [G]  illustrates  the  fact,  noted  earlier, 
that  the  weight  and  c.g.  position  are  defined  separately  by  different  components.  For  example,  ‘x0*  is 
obtainable  from  either  the  normal  force  or  the  side  force  components,  ‘y0*  from  either  side  force  or 
rolling  moment,  and  *zQ*  from  either  normal  force  or  rolling  moment.  Likewise  the  tare  weight  ‘W*  is 
obtainable  from  either  the  normal  force,  side  force,  or  axial  force  components.  There  is,  however,  no 
need  to  evaluate  the  weight  or  tne  c.g.  position  explicitly;  the  objective  of  the  tare  procedure  is  the 
evaluation  of  the  balance  component  load  vector  resulting  from  the  tare  weight,  and  this  is  given  by 
Equation  (53)  in  terms  of  the  matrix  [G]  and  the  vector  made  up  of  trigonometric  functions  of  the 
tare  attitudes. 

4.3  The  ‘Buoyancy’  Condition 

The  hypothetical  condition  of  ‘buoyancy*  can  be  defined  as  that  condition  for  which  the 
electrical  outputs  from  the  various  component  strain  gauge  bridges  result  only  from  resistive  imbalances 
due  to  influences  other  than  applied  loading.  Such  influences  would  include,  for  example,  mismatching 
of  the  individual  strain  gauge  resistance  values,  strains  resulting  from  the  bonding  of  the  strain  gauges 
to  the  flexure  elements,  and  strains  arising  from  assembly  of  multi -part  balances.  The  component  out¬ 
puts  arising  from  any  such  ‘non-load*  influences  will  be  termed  the  ‘buoyant*  component  offsets,  and 
will  be,  by  definition,  independent  of  any  loads  applied  to  the  balance,  as,  for  example,  by  variation  of 
the  balance  attitude  during  the  acquisition  of  tare  data. 

4.4  Determination  of  the  Matrix  [G]  and  the  Buoyant  Offsets 

The  acquisition  of  balance  tare  data  consists  of  determining: 

(a)  the  model/balance  orientation  in  terms  of  the  pitch  (a)  and  roll  (0)  angles  common  to  many 
model  support  systems;  alternatively  equivalent  pitch  and  roll  angles  can  be  established  by 
computing  the  Euler  angles  in  a  ‘yaw-pi tch-roll’  rotation  sequence  which  yield  the  direction 
cosine  matrix  defining  the  orientation  of  the  balance  axes. 

(b)  the  ‘gross*  outputs  from  each  of  the  ‘n’  balance  components  (mT)  where  the  measurements 
are  made  relative  to  the  ‘zero*  of  the  measuring  system. 

Defining  a  vector  (m°)  to  be  the  ‘buoyant’  offsets,  we  can  express  the  net  component  out¬ 
puts  (which  are  proportional  to  the  true  tare  load)  by: 

mN  *  (mT  -  m°)  (56) 


and  this  net  output  vector  is  related  to  the  vector  of  true  component  tare  loads  (FN)  by  Equation  (38). 
Thus: 


mN  =  [Cl]  Fn  +  [C2]  F* 


(57) 


assuming  a  non-linear  balance  calibration. 

Pre-multiplication  of  both  sides  of  this  equation  by  the  diagonal  matrix  composed  of  the 
reciprocals  of  the  diagonal  elements  of  [Cl] ,  reduces  it  to  the  form  of  Equation  (39)  in  which  the 
calibration  coefficient  matrices  are  independent  of  the  measuring  system.  Thus: 

Fn~  -  [XI]  FN  +  [X2]  F* 

FN~  -  [D]  mN 


where 


(58) 
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Pre-multiplication  of  both  sides  of  Equation  (56)  by  [D[  gives: 


[D1  m"  =  [D]  [mr  -  m°] 


pN-  =  pT~  _  po 


and  thus  Equation  (58)  becomes: 


from  which 


Fn~  =  (Ft~  -  F°~)  =  [XI]  Fn  +  (X2)  F* 


Fn  =  [XI]'1  Ft~  -  [XI]-1  F°~  -  [Xir1  [X2]  F* 


Introducing  the  vector  of  ‘buoyant*  offset  loads,  i.e.  the  loads  equivalent  to  the  ‘buoyant’ 
electrical  offsets  (m°)  as  defined  by  a  particular  linear  calibration  matrix,  we  define  it  to  be: 


F°  =  [XI]'1  F° 


and  that  of  the  gross  tare  loads,  which  include  the  ‘buoyant’  offset  loads,  to  be: 


FT  **  [XI]'1  FT" 


then  Equation  (59)  becomes: 


where 


pN  +  po  =  pT  _  [XI]1  [X2]  F* 


F*  =  f(FN) 


Substituting  for  FN  from  Equation  (53)  we  have: 


[sin  a  "1 

cos  a  sin  0  I  +  F°  *  FT 
cos  a  cos  0J 


-  [XI]'1  [X2]  F 


The  acquisition  of  tare  data  at  ‘m’  different  combinations  of  a  and  $  will  result  in  ‘m’ 
equations  of  this  form  which  can  be  expressed  conveniently  in  matrix  notation: 


(01  [T]  +  [F°]  •  [Ft]  -  [XI J'1  [X2]  [F*] 
(n,3)  (3,m)  (n,m)  (n,m)  (n,n)  (n,p)  (p,m) 


(62) 
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Consider  now  the  matrix  [F°]  of  the  ‘buoyant*  offset  loads.  Each  column  is  composed  of 
the  vector  of  offset  loads  (F°)  for  a  particular  tare  attitude  which  is  defined  as: 

F°  =  [XI]'1  [D]  m° 


By  definition  the  vector  of  the  ‘buoyant’  electrical  offsets  (m°)  is  independent  of  the  orienta¬ 
tion  of  the  balance,  and  thus  the  vector  of  offset  loads  (F°)  will  be  independent  of  attitude  also,  pro¬ 
vided  that  the  linear  calibration  coefficients  (i.e.  the  matrices  [XI]  and  [D] )  are  the  same  for  all  tare 
attitudes.  However,  as  has  been  noted  previously,  it  is  quite  common  to  calibrate  a  balance  so  as  to 
define  separate  positive  and  negative  load  coefficients,  and  in  this  case  there  is  no  guarantee  that  the 
matrix  product 

[XI]  1  [D]  =  [Cl]'1 


will  be  identical  for  all  tare  attitudes,  as  these  may  generate  different  distributions  of  the  signs  of  the 
component  loads.  Thus  the  offset  load  vectors  (F°)  for  each  tare  attitude  may  be  different,  despite  the 
fact  that  the  electrical  offset  vectors  (m°)  must  be  identical. 

The  formulation  depends  upon  the  various  F°  being  identical  for  all  tare  attitudes,  a  condi¬ 
tion  which  clearly  is  satisfied  in  the  most  general  case  only  when  the  ‘buoyant*  electrical  offsets  (m°) 
are  zero.  In  practice  this  can  never  be  assured,  and  in  fact  it  is  very  unlikely  to  be  the  case.  However, 
the  assumption  inherent  in  the  formulation  can  be  validated  by  using  an  iterative  computation  proce¬ 
dure  in  which  the  ‘buoyant*  offsets  (m°)  computed  initially  are  subtracted  from  the  tare  data,  which 
are  then  re-processed  to  provide  new  (smaller)  offsets.  This  iteration  can  be  continued  until  the  calcu¬ 
lated  ‘buoyant’  offsets  are  zero,  thus  satisfying  the  assumption  made.  In  this  case,  necessary  only  when 
separate  calibration  coefficients  for  positive  and  negative  load  are  used,  the  actual  ‘buoyant’  offsets 
would  be  given  by  the  summation  of  all  m°  determined  in  the  iterations. 

We  may  therefore  assume  that  the  vector  F°  is  identical  for  all  tare  attitudes;  consequently 
the  matrix  [F°]  is  composed  of  ‘m*  identical  columns.  The  addition  defined  on  the  L.H.S.  of 
Equation  (62)  can  thus  be  accomplished  very  simply  by  matrix  multiplication  of  expanded  [G]  and 
[T]  matrices.  The  [G]  matrix  is  expanded  to  contain  four  columns,  the  last  of  which  is  made  up  of 
the  vector  F°,  and  the  [T]  matrix  is  expanded  to  four  rows  by  addition  of  a  unit  row.  Thus: 

[G]  [T]  +  [F°]  -  [G  [  F°]  £_T  J  -  [Gl]  [Tl] 

(n,3)  (3,m)  (n,m)  (n,4)  (4,m) 


and  Equation  (62)  now  becomes: 

[Gl]  [Tl]  -  [FT]  -  [XI]'1  [X2]  [F*]  (63) 


in  which  [Tl]  and  [FT]  are  known  from  the  tare  measurement!,  and: 

[F‘]  -  f([FN]) 

This  equation  may  be  solved  iteratively  using  as  a  starting  point  the  solution  obtained  by 
neglecting  the  second  order  terms.  Thus  a  first  approximation  of  (Gl  ] ,  or  the  complete  solution  for 
the  case  of  a  linear  balance  calibration,  is  given  ty: 
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(64) 


[Gl]  =  [G  !  F°]  =  [Ft]  [Tl]'1 
(n,4)  (n,3)  (n,l)  (n,m)  (4,m) 


In  order  for  the  matrices  [FT]  and  [Tl]  to  be  conformable  the  number  of  tare  attitudes  ‘m’  must 
clearly  be  equal  to  4,  a  condition  which  also  satisfies  the  requirement  that  [Tl]  be  square  to  permit  its 
inversion.  Thus,  provided  that  the  four  tare  attitudes  are  chosen  to  yield  a  non-singular  matrix  [Tl] , 
the  required  matrix  [Gl]  is  readily  obtainable. 

We  then  have: 


[Gl]  =  [G  !  F°] 


and  the  component  tare  loads  at  attitude  (ap  0*)  are  given  by: 


£N 


[G] 


sin  Oj 

cos  Oj  sin  0j 
cos  cos  0 j 


The  ‘buoyant’  electrical  offsets  (m°)  are  defined  in  terms  of  the  offset  loads  by: 

m°  =  |D|_1  [XI]  F° 


provided  that  the  calibration  did  not  involve  ‘plus-minus’  coefficients.  Otherwise  it  is  necessary  to 
compute  the  m°  vectors  for  each  tare  attitude  using  matrices  [D]  and  [XI]  whose  elements  correspond 
to  the  correct  component  load  sign  distribution  for  each  tare,  and  subtract  these  ‘approximate’  electri¬ 
cal  offsets  from  the  gross  output  vectors  mT.  The  resulting  matrix  of  approximate  ‘net’  outputs  for  all 
tares  is  used  to  compute  new  offset  load  values  which  are  ultimately  reduced  to  zero  by  iteration.  At 
this  point  the  necessary  assumption  that  all  offset  load  vectors  are  zero  is  valid.  The  actual  ‘buoyant’ 
electrical  offsets  are  obtained  as  the  difference  between  the  measured  electrical  outputs  for  one  of  the 
tare  attitudes  and  the  values  of  the  elements  in  the  ‘net’  output  vector,  for  the  same  tare  attitude,  used 
in  the  last  iteration.  The  solution  for  a  linear  calibration  is  then  complete;  for  a  non-linear  calibration 
an  iterative  solution  is  required,  using  this  linear  approximation  as  a  starting  point. 

4.4.1  Iterative  Solution  for  Non-Linear  Balances 

By  neglecting  the  non-linear  terms  initially,  an  approximation  to  the  offset  load  vector  is 

obtained  —  F°0  say,  where  the  subscript  (0,1, - n)  will  be  used  to  indicate  successive  approximations 

in  the  iteration.  Thus  the  matrix  [F°]0,  having  identical  columns  equal  to  F°0,  can  be  generated  and: 


[PN]0  *  [FTlo  -  lF°lo 

where 

[FT]0  -  JX1|- 1  [D]  [mT] 


and  hence  a  first  estimate  of  the  non-linear  force  matrix  ( F*  ]  0  can  be  calculated.  This  gives  a  further 
approximation,  [FT] ,  say,  to  the  R.H.8.  of  Equation  (63): 


-31  - 


[FT],  =  [Ft]0  -  [XI]1  [X2]  [F*]0 


and,  from  Equation  (64),  we  then  have: 

[Gl] i  =  [FT],  [Tl]-‘ 

The  last  column  of  [Gl]  yields  a  further  approximation  of  the  offset  load  vector  F°!,  and 
hence  the  matrix  [F°]  1?  and  the  iteration  proceeds. 

In  general: 

[FT]n  =  [FT]0  -  [XI]'1  [X2]  [F*](n-i) 

where 

[F*](n-D  =  fllF1]^.!)  -  [F0](n_i)}  =  f{[FN](n-i)} 

[Gl]n  =  [FT]n  [Tl]1 

and  t F° ]  ( n - 1)  is  composed  of  identical  columns  equal  to  the  last  column  of  [Gl ]  ( n > i ) . 

Convergence  of  the  solution  is  established  when: 

[FT]n  -  [FT](n_n  <  5 


where  the  vector  6  consists  of  specified  limits  for  each  balance  component.  The  solution  has  been 
found  to  converge  rapidly  for  typical  non-linear  balances,  the  convergence  being  aided  considerably  by 
the  fact  that  the  loads  involved  in  tare  measurements  are  normally  small  relative  to  the  balance  rated 
loads.  Thus  non-linear  effects  will  be  quite  small  even  when  the  non-linear  calibration  coefficients  are 
relatively  large.  As  a  precaution  against  divergence  of  the  solution,  the  values  of  [FT]n  may  be  com¬ 
pared  with  test  values  related  to  the  component  load  capacities. 

As  in  the  case  of  the  linear  solution,  when  ‘plus-minus’  balance  calibration  coefficients  are 
employed  it  is  necessary,  upon  convergence  of  the  non-linear  solution,  to  iterate  the  entire  procedure 
in  order  to  reduce  the  offset  load  vector  (F°)  to  zero. 

4.6  Application  of  the  Tare  Procedure 

The  objective  of  the  tare  computation  is  the  evaluation  of  the  matrix  [G]  in  Equation  (53) 
and  the  vector  of  ‘buoyant’  electrical  offsets  m°,  both  of  which  are  defined  by  the  matrix  [Gl]  deter¬ 
mined  by  the  procedure  given.  By  using  the  ‘buoyant’  offsets  (m°),  which  are  defined  in  the  measure¬ 
ment  units  of  the  data  acquisition  system,  all  loads  are  automatically  computed  relative  to  absolute 
zero  load,  and  thus  support  rotations  and  the  deflected  orientation  of  the  balance  axes  system  are 
readily  obtainable.  Determination  of  the  equivalent  pitch  and  roll  angles,  (as  the  second  and  third 
Euler  rotations  in  a  *yaw-pitch-roir  sequence  of  rotations  which  yields  the  direction  cosine  matrix  of 
the  required  axes  orientation),  then  enables  direct  use  of  Equation  (63)  to  establish  the  balance  com¬ 
ponent  loads  of  the  tare  weight.  These  are  subtracted  from  the  total  component  loads  to  arrive  at  the 
required  aerodynamic  loading. 
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It  is  clearly  unnecessary  to  evaluate  the  tare  weight  and  centre  of  gravity  position  explicitly. 
However,  as  indicated  in  Section  4.2,  the  elements  of  the  matrix  [G]  can  be  defined  analytically  in 
terms  of  W,  xQ,  y0,  and  zQ  for  any  particular  balance,  from  which  it  is  possible  to  evaluate  these  quan¬ 
tities  by  substitution  of  the  numeric  values  of  the  matrix  elements.  To  do  so  often  provides  a  useful 
‘diagnostic’  indication  of  the  validity  of  the  data;  obviously  incorrect  results  usually  indicate  a  violation 
of  the  assumption  that  only  gravitational  forces  are  acting  on  the  balance,  perhaps  as  a  result  of  ‘fouling’ 
between  the  ‘live’  and  ‘ground’  sides  of  the  model/balance  system. 

In  situations  combining  a  relatively  heavy  model  and  a  relatively  flexible  support,  it  may  be 
desirable  to  correct  the  tare  attitudes  for  rotations  of  the  support  under  load.  Usually  the  support 
stiffness  constants  will  be  defined  in  terms  of  the  balance  component  loads,  and  thus  the  rotations  are 
easily  computed  from  the  balance  components  of  the  tare  load.  Following  application  of  the  tare  pro¬ 
cedure  to  evaluate  the  matrix  [G] ,  the  balance  tare  components  can  be  determined  from  Equation  (53), 
the  support  rotations  calculated,  and  the  tare  attitudes  corrected  accordingly.  The  tare  procedure  can 
then  be  applied  again,  using  the  corrected  attitudes,  and  [G]  and  m°  obtained.  Such  a  correction  for 
support  rotations  is  essentially  exact,  as  the  change  in  the  computed  component  tare  loads  as  a  result 
of  the  deflection  will  be  very  small  compared  to  the  magnitude  of  the  component  loads. 

4.6  Optimization  of  the  Tare  Attitudes 

As  described,  the  tare  procedure  requires  data  to  be  recorded  at  four  unique  model/balance 
attitudes,  though  note  here  that  it  is  the  orientation  of  the  balance  axes  system  rather  than  the  model 
which  is  important.  This  represents  the  most  general  case  where  the  balance  system  is  responsive  to  all 
three  components  of  the  tare  weight,  (Wx,  Wy,  and  Wz);  abbreviations  in  which  one  or  more  of  these 
components  do  not  influence  the  balance  readings  will  be  discussed  later.  The  question  which  now 
arises  is,  what  combination  of  the  four  tare  attitudes  will  provide  the  most  accurate  result? 

A  possible  answer  is  suggested  by  observation  of  Equation  (64)  which  states: 

[Gl]  =  [Ft]  [Tl]-1 


It  is  clear  that  the  combination  of  attitudes  must  result  in  a  non-singular  matrix  [Tl] ,  which  consists 
of  trigonometric  functions  of  the  pitch  and  roll  angles  defining  the  balance  orientation.  Further,  as  the 
matrix  [FT]  is  composed  of  measured  balance  data  which  are  subject  to  experimental  error,  minimizing 
the  absolute  values  of  the  elements  of  the  inverse  of  [Tl]  might  logically  be  expected  to  minimize  the 
effect  of  these  experimental  errors  following  the  multiplica  lion  indicated.  This  condition  would  be 
equivalent  to  maximizing  the  absolute  value  of  the  determinant  of  the  matrix  [Tl] . 

Using  this  criterion,  the  optimum  combination  of  tare  attitudes, 


i  =  1,4] 


were  determined  using  an  existing  computer  program  based  on  the  Simplex  method1 The  result 
obtained  is  summarized  below: 


Tare  #  (i) 

“i 

0i 

1 

ai 

01 

2 

(0,  +  180°) 

3 

<*2 

(0,  +  90°) 

4 

<*2 

(01  +  270°) 
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The  value  of  <px  is  arbitrary,  the  maximum  absolute  value  for  the  determinant  occurring  when: 

cq  =  “a2  =  sin’1  (l/y/3)  —  35.3° 


and 

|A|  -  16/(3\/3)  -  3.08 


Using  this  combination  of  attitudes,  the  analytic  expression  for  the  determinant  simplifies  to: 


A  =  4  cos  cq  cos  a2  (sin  oq  -  sina2) 


Figure  2  illustrates  the  variation  of  |A|  with  the  total  range  of  the  pitch  angle  (a),  and  the  mid  point  of 
this  range. 


It  is  suggested  that  the  arbitrary  datum  roll  angle  (0j )  should  be  chosen  such  that  all  balance 
components  influenced  by  changes  of  roll  angle  should  experience  finite  loads  at  all  tare  attitudes.  This 
will  ensure  a  unique  choice  of  ‘plus-minus’  calibration  coefficients  for  the  component  load  sign  distri¬ 
bution  of  each  tare  attitude.  A  datum  roll  angle  of  45  degrees,  combined  with  as  large  a  range  of  the 
pitch  angle  as  practicable  (up  to  the  optimum  value  indicated  above),  would  therefore  be  appropriate. 

4.7  Abbreviations  of  the  General  Case 

Situations  can  arise  in  which  one  or  more  of  the  balance-axes  tare  weight  components,  (Wx, 
Wy,  and  Wz),  may  produce  no  variation  of  any  balance  component  load  as  the  model/balance  attitude 
is  changed.  Such  would  be  the  case,  for  example,  with  a  balance  designed  to  measure  axial  force  only,  as 
here  the  components  Wy  and  Wz  have  no  influence  on  the  balance  output.  In  such  cases  Equation  (53) 
can  be  abbreviated  by  elimination  of  the  appropriate  balance-axes  tare  load  component,  or  components. 

Another  situation  requiring  abbreviation  of  the  general  formulation  arises  when  the  model 
attitude  variation  is  restricted  to  either  pitch  or  roll  for  some  practical  consideration.  For  example, 
should  there  be  no  capability  to  alter  the  roll  angle,  it  would  not  be  possible  to  acquire  four  independ¬ 
ent  tares  and  the  result  would  be  a  singular  matrix  [Tl] ;  here  the  elements  of  the  second  and  third 
rows,  (i.e.  for  Wy  and  Wz),  would  not  be  independent.  To  obtain  a  non-singular  matrix  [Tl]  in  this 
case,  it  is  necessary  to  eliminate  either  Wy  or  Wz  from  the  formulation;  each  balance-axes  tare  compo¬ 
nent  eliminated  reduces  the  required  number  of  tares  by  one. 

There  are  consequently  two  situations  requiring  the  use  of  an  abbreviated  form  of  the 
formulation.  In  the  first  case  described,  there  is  no  loss  of  generality  as  the  balance-axes  components 
are  eliminated  only  because  they  do  not  affect  the  balance  measurements.  However,  in  the  second  case 
the  component  must  be  eliminated  as  it  is  not  independent.  For  whichever  reason,  elimination  of  a 
component  brings  about  a  contraction  of  the  [G]  matrix,  the  contracted  form  correctly  relating  the 
balance  tare  components  and  the  balance-axes  (gravity)  components  of  the  tare  load  only  when  the 
attitude  is  compatible  with  the  tare  attitude.  Thus,  if  the  tare  data  are  recorded  at  constant  roll  angle, 
the  contracted  [G]  matrix  calculated  by  the  tare  procedure  will  give  correct  balance  tare  components 
only  at  that  constant  roll  attitude. 

The  most  general  tare  condition  is  that  given  in  Equation  (52);  the  six  possible  abbreviated 
forms  are  noted  below: 

(a)  Wx  eliminated:  FN  =  W[G2  cos  ot  sin  0  +  G3  cos  a  cos  0] 

(b)  Wy  eliminated:  FN  =  W[Glsina  +  G3  cos  a  cos  <*>] 
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(c) 

Wz  eliminated: 

Fn  =  W[Glsina  +  G2  cos  a  sin  01 

(d) 

Wy  and  Wz  eliminated: 

Fn  =  W[G1  sin  a] 

(e) 

Wx  and  Wz  eliminated: 

Fn  =  W[G2  cos  a  sin  0] 

(f) 

Wx  and  Wy  eliminated: 

Fn  =  W[G3  cos  a  cos  0] 

Forms  (a),  (b),  and  (c)  require  three  tares  while  (d),  (e),  and  (f)  require  two  tares  only. 

Generally  the  latter  three  abbreviations  would  be  used  when  the  balance  was  influenced  by 
one  gravity  component  only,  and  thus  are  particular  to  balances  designed  to  measure  loads  in  just  one 
plane.  Cases  (a),  (b),  and  (c)  represent  the  more  interesting  situation  brought  about  by  physical  restric¬ 
tions  in  the  model  attitude  control. 

Case  (a)  corresponds  to  all  tares  being  acquired  at  constant  pitch  angle.  The  contracted  [G] 
matrix  contains  two  columns,  appropriate  to  the  gravity  components  Wy  and  Wz,  which  correctly 
define  the  balance  tare  components  at  the  given  constant  pitch  angle  only.  However,  the  computed 
electrical  offsets  (m°)  will  be  ‘buoyant’  only  if  the  constant  pitch  angle  is  0°.  At  any  other  constant 
value  the  offset  loads  calculated  for  each  balance  component  will  be  the  sum  of  the  true  ‘buoyant’ 
offset  load  and  the  component  load  generated  by  the  ‘Wx’  component  of  the  tare  load.  The  offset  load 
for  rolling  moment  will  consequently  be  ‘buoyant’  regardless  of  the  value  of  the  pitch  angle,  and, 
knowing  the  tare  weight  and  centre  of  gravity  location,  the  true  ‘buoyant’  offsets  for  the  other  compo¬ 
nents  could  be  determined  if  desired.  A  similar  situation  exists  for  the  abbreviated  cases  (e)  and  (f), 
while  for  case  (d)  the  calculated  offset  for  the  ‘X’  balance  component  is  ‘buoyant’. 

Cases  (b)  and  (c)  both  correspond  to  the  acquisition  of  all  tares  at  constant  roll  angle.  A  roll 
angle  of  0  or  180  degrees  requires  the  use  of  case  (b),  while  case  (c)  is  appropriate  if  the  angle  is  90  or 
270  degrees.  For  any  angle  other  than  those  noted  either  abbreviation  may  be  used,  although  the 
logical  choice  would  be  to  eliminate  the  gravity  component  exhibiting  the  smaller  variation  with  pitch 
angle.  For  these  cases  the  offsets  computed  are  the  true  ‘buoyant’  offsets,  regardless  of  the  value  of 
the  constant  roll  angle. 

The  optimum  tare  attitude  combinations  for  these  abbreviated  cases,  again  based  on  maxi¬ 
mizing  the  absolute  value  of  the  determinant  of  [Tl] ,  are  given  below: 


Case 

a2 

±L 

02 

03 

|A|  max  at 

(a) 

0° 

0° 

0° 

01 

0° 

-01 

0,  =  120° 

(b) 

“l 

0° 

-«1 

0°/180° 

0°/180° 

0°/180° 

a,  =  120° 

(c) 

“l 

0° 

-«1 

90°/270° 

90°/270° 

90°/270° 

a,  =  120° 

(d) 

«1 

-“1 

01 

01 

— 

a,  =  90° 

(e) 

0° 

0° 

01 

-01 

0,  =  90° 

(f) 

0° 

0° 

01 

0,  + 180° 

— 

0,  =  0° 

5.0  CONCLUDING  REMARKS 

This  report  has  discussed  the  calibration  and  application  of  strain  gauge  balances  with  par¬ 
ticular  reference  to  their  use  in  wind  tunnels.  The  procedures  for  multi-component  loading  calibrations 
and  the  application  of  force  balance  tares  utilizing  the  concept  of  ‘buoyancy’,  were  developed  at  NAE 
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as  part  of  a  program  aimed  at  rationalizing  the  balance  calibration  and  data  reduction  procedures  used 
in  connection  with  tests  performed  in  the  5  ft.  X  5  ft.  Wind  Tunnel  and  the  15  in.  X  60  in.  Two- 
Dimensioned  Test  Section.  The  procedures  described  have  been  in  use  at  NAE  for  a  considerable  period 
and  have  proven  to  be  very  satisfactory.  Procedures  for  solution  of  the  second  order  balance  equations 
have  been  included  to  make  this  report  a  fairly  complete  description  of  methods  for  calibration  and 
use  of  strain  gauge  balances.  However,  the  non-iterative  method  outlined  is  not  employed  at  NAE  for 
reduction  of  balance  data  from  the  5  ft.  X  5  ft.  Wind  Tunnel. 

A  general  method  for  the  calibration  of  multi-component  strain  gauge  balances  has  been 
presented,  and  compared  with  the  common  procedure  in  which  the  components  are  loaded  independ¬ 
ently  and  constant  secondary  loads  are  used  to  determine  the  non-linear  interactions.  In  presenting 
both  methods  the  equations  defining  the  balance  behaviour  have  been  limited  to  second  order  terms 
as,  from  the  experience  of  many  users,  higher  order  effects  have  been  shown  to  be  negligible  for  typical 
balances.  However,  both  methods  are  easily  extendable  to  include  higher  order  effects  if  desired. 

A  major  difference  between  the  two  methods  is  that  the  latter  specifically  limits  the  order  of 
the  interaction  terms,  usually  to  second  order,  by  restricting  the  number  of  component  loads  acting 
simultaneously,  whilst  the  former  allows  the  existence  of  higher  order  effects  but  limits  the  terms  in 
the  balance  equations  by  the  choice  of  the  order  of  the  polynomial  to  which  the  data  are  fitted.  Any 
terms  in  the  general  balance  equation  which  are  deemed  to  be  negligible  in  a  particular  situation  may 
be  eliminated,  (the  interaction  of  a  small  axial  force  on  a  large  normal  force  for  example),  thus  abbre¬ 
viating  the  equations  and  hence  the  number  of  loading  cases  necessary  to  obtain  a  solution.  Note  how¬ 
ever  that,  if  the  component  loads  involved  remain  non-zero  during  the  calibration  process,  (normal  and 
axial  forces  in  the  example),  elimination  of  the  term  from  the  equation  does  not  remove  the  actual 
interaction,  and  the  abbreviated  formulation  still  represents  a  best  solution  of  the  equations  when 
using  the  actual  data. 

The  proposed  method  allows  standardization  of  calibration  procedures  for  both  ‘internal’ 
and  ‘external’  types  of  balance.  The  use  of  a  single  varying  calibration  load  permits  the  use  of  loadings 
which  are  representative  of  those  to  be  expected  in  a  test  environment,  thus  ensuring  reasonable 
similarity  of  the  deflections  occurring  in  both  calibration  and  use  of  the  balance.  All  calibration  co¬ 
efficients  can  be  obtained  without  the  need  of  ‘secondary’  loads,  thus  simplifying  the  calibration 
set-up.  The  formulation  does,  however,  permit  the  inclusion  of  constant  ‘secondary’  loads  on  one  or 
more  components,  thus  providing  a  means  of  accounting  for  the  weight  of  the  calibration  body  so 
that  the  load  on  each  component  may  be  specified  relative  to  absolute  zero  load  —  i.e.  a  ‘buoyant’ 
zero.  In  combination  with  the  method  for  defining  ‘buoyant’  tares,  this  avoids  ‘initial  load  effects’  in 
the  use  of  non-linear  balances. 

The  choice  of  a  ‘linear’  versus  a  ‘non-linear’  calibration  should  be  based  upon  whether  or  not 
the  simpler  linear  matrix  can  provide  adequate  accuracy  for  the  particular  test  requirements.  With  the 
multi-component  loading  technique  a  linear  solution  is  really  a  linearization  of  the  complete  balance 
behaviour.  This  may  provide  a  better  accuracy  than  the  linear  solution  obtained  using  the  technique  of 
loading  the  components  independently,  as  in  this  case  no  non-linear  effects  are  present  during  the 
calibration.  When  ‘linearizing’  a  non-linear  system  it  should  be  noted  that  the  linear  curve  fit  of  the 
data  should  be  forced  to  pass  through  the  actual  data  value  at  zero  load;  performing  the  curve  fit  with 
a  ‘free’  intercept  will  result  in  a  uni-directional  bias  upon  subtraction  of  ‘tares’. 

The  ability  of  the  derived  matrix  of  sensitivity  and  interaction  coefficients  to  adequately 
represent  the  actual  balance  behaviour  should  be  verified,  by  using  it  to  re-calculate  the  applied  load 
for  some  check  cases  which  closely  simulate  the  expected  test  loadings.  A  convenient  means  of  express¬ 
ing  the  result  is  to  compute  the  percentage  error  in  the  magnitude  of  the  load  vector,  and  the  angular 
error  in  its  direction. 

The  most  satisfactory  representation  of  the  balance  behaviour  will  probably  be  obtained  by 
treating  the  response  to  positive  and  negative  loads  separately.  It  is,  however,  most  important  to 
recognize  that  it  is  the  sign  of  the  strain  in  the  flexure  element  which  matters,  not  the  sign  of  the  load 
for  some  ‘derived’  component  which  is  a  function  of  two  or  more  flexure  loads.  It  seems  superfluous 
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to  state  this,  but  the  author  is  aware  of  a  case  in  which  the  behaviour  of  a  6-component  balance  was 
specified  by  all  linear,  load  squared,  second  order  load  cross-product,  load  cubed,  and  load  ‘absolute 
value'  terms;  however  the  absolute  value  terms  were  based  on  the  signs  of  derived  component  loads,  a 
fact  which  considerably  undermined  any  confidence  which  might  have  been  placed  in  the  inclusion  of 
the  load  cubed  terms! 

A  non-iterative  method  for  solution  of  the  second  order  balance  equations  has  been  pre¬ 
sented.  However,  it  is  felt  that  the  approximation  involved  cannot  be  justified,  as  the  iterative  solution 
converges  rapidly  for  typical  balances.  Moreover,  the  non-iterative  method  cannot  accommodate  the 
use  of  separate  calibration  coefficients  appropriate  to  positive  and  negative  component  loads.  Two 
methods  of  accounting  for  such  ‘plus/minus'  calibration  matrices  have  been  described.  One  method 
involves  the  synthesis  of  a  specific  matrix,  appropriate  to  a  given  vector  of  component  loads,  from  the 
general  ‘plus/minus’  matrix.  This  specific  matrix  must  be  inverted  as  part  of  the  solution  for  the  com¬ 
ponent  loads.  The  second  method  uses  pre-defined  coefficient  matrices  relating  to  the  ‘signed’  and 
‘absolute’  values  of  the  component  loads,  and  in  this  case  no  matrix  inversion  is  required  although 
additional  matrix  multiplications  are  necessary.  From  a  computational  efficiency  standpoint,  the 
latter  method  is  to  be  preferred. 

The  procedure  described  for  the  application  of  balance  tares  uses  a  concept  of  ‘buoyancy’  to 
overcome  the  iterative  nature  of  the  problem;  this  arises  if  the  attitude  is  not  sensed  within  the  model 
and  must  therefore  be  corrected  for  support  deflections  under  load.  Additional  advantages  of  the 
method  are  the  determination  of  the  ‘buoyant’  electrical  offsets  for  each  component,  which  provides 
a  convenient  means  of  monitoring  the  long-term  stability  of  the  balance,  and  the  fact  that  the  derived 
weight  and  centre  of  gravity  information  can  provide  a  useful  ‘diagnostic’  of  the  operation  of  the  bal¬ 
ance  —  e.g.  freedom  from  ‘fouling’  between  the  ‘live’  and  ‘ground’  members.  The  method  requires 
measurement  of  the  component  outputs  and  the  balance  orientation  at  a  maximum  of  four  independ¬ 
ent  positions.  The  orientation  is  defined  in  terms  of  the  pitch  and  roll  angles  available  with  a  common 
type  of  model  support  system,  and  optimum  combinations  of  these  angles  are  suggested.  For  support 
systems  which  do  not  utilize  the  specified  pitch  and  roll  angles  for  model  positioning,  the  necessary 
angles  can  be  determined  from  the  matrix  of  direction  cosines  which  defines  the  orientation  of  the 
balance  axes  system. 

The  use  of  the  ‘buoyant’  tare  approach,  in  combination  with  a  definition  of  loads  on  an 
absolute  scale  during  balance  calibration,  avoids  the  problem  of  initial  load  effects  in  the  use  of  non¬ 
linear  balances. 
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FIG.  1:  SIGN  CONVENTIONS  FOR  FORCES  AND  DISTANCES 


ABSOLUTE  VALUE  OF  THE  DETERMINANT  OF  [T]  (|AI) 
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FIG.  2:  VARIATION  OF  THE  ABSOLUTE  VALUE  OF  THE  DETERMINANT  OF  [T] 
WITH  PITCH  ANGLE  (cr>  FOR  THE  OPTIMIZED  ATTITUDES 
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